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Abstract 

The construction of the effective Lagrangian relevant for the mesonic 
sector of QCD in the large limit meets with a few rather subtle problems. 
We thoroughly examine these and show that, if the variables of the effec- 
tive theory are chosen suitably, the known large Nc counting rules of QCD 
can unambiguously be translated into corresponding counting rules for the 
effective coupling constants. As an application, we demonstrate that the 
Kaplan-Manohar transformation is in conflict with these rules and is sup- 
pressed to all orders in The anomalous dimension of the axial singlet 
current generates an additional complication: The corresponding external 
field undergoes nonmultiplicative renormalization. As a consequence, the 
Wess-Zumino-Witten term, which accounts for the 11(3)^x11(3)^ anoma- 
lies in the framework of the effective theory, contains pieces that depend 
on the running scale of QCD. The effect only shows up at nonleading order 
in but requires specific unnatural parity contributions in the effective 
Lagrangian that restore renormalization group invariance. 



Pacs: 11.15.Pg, ll.30.Rd, 12.38.Aw, 12.39.Fe, 14.40.-n 

Keywords: QCD, 1/Nc expansion, effective Lagrangian, chiral symmetry, 
anomalies, chiral perturbation theory 



Work supported in part by Schweizerischer Nationalfonds 



Contents 



1 Introduction 



2 Effective action 



3 QCD at large 



4 Relation between QCD and gluodynamics 

5 Massless quarks 



6 Simultaneous expansion in p and 



7 Effective theory in the U(l) sector 

8 Extension to U(3l 



9 Chiral symmetry 



10 Full effective Lagrangian 



11 Potential^ 



12 Dependence on the vacuum angle 



13 Kaplan-Manohar transformation at large 



14 Renormalization of the operators A" and u: 



15 Renormalization of the effective action 



3 
4 
5 
7 
9 
10 
12 
14 
15 
17 
19 
21 
23 
24 
26 

28 
31 
33 
36 
38 
39 
40 



16 Dependence of the effective coupling constants on the running 
scale of QCD 



17 Higher orders and loops 



18 Matching U(3) and SU(3) 



19 Anomalies 



20 Wess-Zumino-Witten term 



21 Unnatural parity part beyond leading order 



22 Radiative transitions 



1 



23 Summary and conclusion 



A Construction of the effective Lagrangian 



B Full second order SU(3) Lagrangian 



C Renormalization of the WZW term 



41 
44 
50 
53 



2 



1 Introduction 



The low energy properties of QCD are governed by an approximate, sponta- 
neously broken symmetry, which originates in the fact that three of the quarks 
happen to be light. If m„,mrf,ms are turned off, the symmetry becomes exact. 
The spectrum of the theory then contains eight massless pseudoscalar mesons, 
the Goldstone bosons connected with the spontaneous symmetry breakdown. 

If the number of colours is taken large, the quark loop graph which gives rise 
to the U(l)-anomaly is suppressed This implies that, in the limit — > oo, 
the singlet axial current is also conserved: The theory in effect acquires a higher 
degree of symmetry. Since the operator qq fails to be invariant under the extra 
U(l)-symmetry, the formation of a quark condensate, {0\qq |0) 7^ 0, implies that 
this symmetry is also spontaneously broken 0. The spectrum of QCD, therefore, 
contains a ninth state, the rj', which becomes massless if not only mu,md,ms are 
turned off, but if in addition the number of colours is sent to infinity. 

Chiral symmetry imposes strong constraints on the properties of the Gold- 
stone bosons. These may be worked out in a systematic manner by means of 
the effective Lagrangian method, which describes the low energy structure of 
the theory in terms of an expansion in powers of energies, momenta and quark 
masses |]. The fact that, in the large limit, the rj' also plays the role of a 
Goldstone boson implies that the properties of this particle are subject to analo- 
gous constraints, which may again be worked out by means of a suitable effective 
Lagrangian. 

The main features of the effective theory relevant for the meson sector of QCD 
in the large limit were discovered long ago [Q. The systematic analysis in the 
framework of chiral perturbation theory was taken up in ref. [Q, where the Green 
functions of QCD were studied by means of a simultaneous expansion in powers 
of momenta, quark masses and A considerable amount of work concerning 

the structure of QCD at large A^ has been carried out since then and 
quite a few phenomeno logical applications of the 1/A^ expansion have appeared 
in the literature. For a review of these, in particular also for a discussion of the 
rj-rj' mixing pattern, we refer to [^. Moreover, the large A^ properties of the 
effective theory relevant for the baryons were recently investigated in detail |]23| . 

In the present paper, we examine the foundations of the effective theory in 
the sector with baryon number B = and show that this leads to new insights 
into the low energy structure of QCD in the large A^ limit. The framework is 
more complicated than in the case of A^ = 3, because an additional low energy 
scale appears, related to the mass of the 77'. In order to firmly establish our 
claims, we first need to put the effective theory on a solid basis and demonstrate 
that the effective coupling constants necessarily obey the counting rules indicated 
in ref. (sections @-12 and appendix ^). Using these, we then show that the 
transformation introduced by Kaplan and Manohar is in conflict with the large 
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properties of QCD: The parameter A occurring in the transformation 



(cychc 



u 



d 



u) 



'l.V. 



vanishes to all orders of the expansion (section |T3|). Next, we discuss the con- 
sequences of the fact that the dimension of the singlet axial current is anomalous 
and determine the corresponding dependence of the effective coupling constants 
on the running scale of QCD (sections |l^l6). The matching between the ef- 
fective theories relevant for a finite and an infinite number of colours is worked 
out in section ^ and appendix |B|. Finally, in sections [T9H2I1 and appendix 0, 
we examine the modifications required to extend the Wess-Zumino-Witten term 
from SU(3)«xSU(3)i to U(3)«xU(3)i. 

Some of the results described in the present paper were reported earlier 
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The application to the masses, decay constants and photonic transitions of the 
pseudoscalar mesons will be given elsewhere E^. A further application concerns 



the low energy properties of the correlation function 



X{q^) = -I dxe"i%Q\Tuo{x)uo{Q) |0) 



167r2 



1.21 



where we have absorbed the coupling constant in the gluon field. In particular, 
the manner in which the topological susceptibility x(0) and the derivative x'(0) 
depend on the light quark masses is quite remarkable [E^ p5 



2 Effective action 

Our analysis is based on the known large behaviour of QCD. We work with 
the effective action of this theory, which describes the response of the system to 
the perturbation generated by a set of external fields, 

^QCD = >C°cD +^7^(^M + -^{s - i-f5p)q -Ouj . (2.1) 

The term C^^d describes the limit where the masses of the three light quarks 
and the vacuum angle are set to zero. The external fields f^(x), a^{x), s{x), p{x) 
represent hermitean 3x3 matrices in flavour space. The mass matrix of the three 
light quarks is contained in the scalar external field s{x). The vacuum angle d{x) 
represents the variable conjugate to the operator a;(x) specified in eq. (pT^). In 
Euclidean space, the integral 




is the winding number of the gluon field, so that uj{x) may be viewed as the 
winding number density. 
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The effective action represents the logarithm of the vacuum-to- vacuum tran- 
sition amphtude, 

gi5e^{w,P,e} ^ /Qo^^lOin) (2.2) 

\ I ' v,a,s,p,a \ / 

The coefficients of the expansion of Seff{v,a, s,p,6} in the external fields are 
the connected correlation functions of the vector, axial, scalar and pseudoscalar 
quark currents and of the winding number density, in the massless theory. Lorentz 
invariance implies that these can be decomposed into scalar functions, with coef- 
ficients that contain the external momenta and the tensors g^i,, e^i^p„. In view of 
the fact that the square of e^j^po- can be expressed in terms of Qfj^i^, there are two 
categories of contributions: The natural parity part of the effective action, which 
collects the pieces that do not contain the e-tensor, and the unnatural parity part, 
where this tensor occurs exactly once. 

The consequences of the symmetry properties of QCD for the natural parity 
part of the effective action are remarkably simple: The Ward identities related 
to the conservation of the vector and axial currents imply that this part of the 
effective action is invariant under chiral gauge transformations, 

< = V,r,V^ + tVnd.V^ , /; = VJ,V^ + tV,d,V^ (2.3) 
s' + ip' = Vn{s + ip)V^ , e' = 6' + ilndetVR -ilndetV^ , 

with T/, = Vf, + a^, l^ = Vf^- and Vii{x), Vl{x) e U(3). 

Chiral U(l) transformations play an essential role for the properties of the 
theory at large A^. It is important that our list of external fields includes a 
source for the singlet axial current: the trace of the matrix a^. The divergence of 
this current contains an anomaly proportional to lj. The above transformation 
law automatically accounts for this term, through the shift in the vacuum angle 
that is generated by a U(l) rotation. 

The U(l) anomaly is not the only one relevant in our context, but the remain- 
ing anomalies of the effective action are independent of the interaction, so that 
they only affect the unnatural parity part. We first investigate the low energy 
structure of the natural parity part, which, as stated above, is gauge invariant: 
Seff{v', a', s',p', 9'}np — Seff{v, a, s,p, 9}„p. The unnatural parity part - in partic- 
ular, the Wess-Zumino-Witten term - only contributes at nonleading orders of 
the low energy expansion and will be discussed in detail later on. 



3 QCD at large Nc 

The well-known leading logarithmic formula 

7J^ = R 1 ( 2/A2 ^ ' (ll^c - 27V^) 

(47r)2 /3oln(/xVA2^j3) 3 
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implies that the running couphng constant tends to zero when Nc becomes large, 
~ 1/^- At leading order of the expansion, the Green functions are 

dominated by those graphs that contain the smallest possible number of quark 

loops. For the correlation functions of the quark currents, graphs with one such 

loop generate the leading contributions. 

Consider the connected correlation function formed with Uj quark currents 

ji — q FiQ and riui winding number densities 

Gn,n^ = (0| Tji{xi) ■ ■ ■jn^{xn^)uj{yi) ■ ' ' uj{ynj \0) c (3-2) 

and denote the fraction due to graphs with i quark loops by G^.^^. The large 
counting rules of perturbation theory imply that this quantity represents a term 
of order 

<.n. = 0(Ar2-^--^), £^0,1,... (3.3) 

The leading power is independent of the number rij of quark currents, but de- 
creases with n^. This implies that the dependence of the effective action on the 
vacuum angle is suppressed. The leading terms in the expansion of the 
effective action may be characterized by the formula 

Seff ^ N^^ So{i9} + N,Si{v,a,s,p,i9} + S2{v,a,s,p,i9} + . . . , (3.4) 

where Se collects the contributions from graphs with i quark loops and 

(3.5) 

Since the external fields v, a, s,p contained in 5*1 are attached to one and the same 
quark loop, the expansion of Si in terms of these fields generates expressions that 
involve a single trace over the fiavour indices. Likewise, 5*2 represents a sum of 
contributions containing at most two fiavour traces, etc. 

The dominating contributions to the effective action arise from those graphs 
that do not involve quark lines. Their sum represents the effective action of 
gluodynamics, which only depends on the vacuum angle. The expansion in powers 
of 9 yields the connected correlation functions of the field u;{x), 

S^f] = ^Jdxidx2e{xi)e{x2){0\Tu;{xi)u;{x2)\0)au 

~i[ J^^^ ' ' •'^^4 6'(a;i) • • ■ 9{x4){0\Tu!{xi) ■ ■ ■u!{x4)\0)ao + • • ■ 

The leading contribution in the l/Nf. expansion of the effective action of gluody- 
namics coincides with the one occurring in QCD, 

S^f]^N,'Som + 0{N,). (3.6) 



6 



In particular, the topological susceptibility of gluodynamicsQ 



= ^^° = jdx{0\Tuj{x)uj{0)\0) 



GD ; 



represents a term of order 1. 

Since gluodynamics possesses a mass gap (excitation energy of the lightest 
glueball), the correlation functions of u decay on a scale of order Aqcd- We are 
analyzing the theory for external fields that vary slowly on this scale. In that 
regime, we may expand the factor 9{xi) . . . 9{xn) around the point xi and express 
the functional through a derivative expansion: 

So{^} = Jdx{-eoi^) + d^^d''^eii^) + ...} , (3.7) 

where eo{'&) , ei{^) , . . . are ordinary functions of a single variable. The quantity 
N^eo{^) represents the energy density of the vacuum for the case where, in the 
vicinity of the point under consideration, the vacuum angle is taken constant. 
The term N^d^d d^d ei{d) describes the change in the energy density that arises 
if the vacuum angle has a nonvanishing gradient there, etc. 



4 Relation between QCD and gluodynamics 

The large counting rules imply that the quantities eo('(9), eii^), . . . , which 
occur in the derivative expansion ( |3.7D of the effective action of gluodynamics, 
depend on the variables 9 and only through the ratio d = 9/Nc. This is 
puzzling, because the effective action of QCD is periodic in 9 with period 27i. 
The resolution of the paradox is given in ref. ||2^, where the relation between 



QCD and gluodynamics is discussed in detail. The paradox is related to the fact 
that the partition function Zq^ involves a sum over all gluon field configurations, 
while Zqcu only extends over those on which fermions can live. In gluodynamics, 
the sum over all configurations involves fractional winding numbers - only the 
quantity Nc v must be an integer. In the presence of fermions, however, the 
boundary conditions imply that v itself must be an integer. Even if their mass is 
sent to infinity, the fermions thus exert a restriction on the gluon configurations 
to be summed over: Only configurations with integer winding number contribute. 
This implies that, if either the number of colours or the quark masses are sent 
to infinity, the QCD partition function does not in general tend to the partition 



^As written, the formula refers to Euclidean space. In Minkowski space, the expression reads 
Tgd = -i\dxif)\Tuj{x)uj{^)\G)G,, . 
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function of gluodynamics. Instead, it approaches the one obtained by restricting 
the sum over all gluon configurations to those with integer winding number, given 
by the projection 

k=0 

which removes fractional winding numbers. While Zqq{6) is periodic in 6 with pe- 
riod 27iNc, the projection Zqu{9) has the 27r periodicity characteristic of Zqcd(^)- 
If 9 is taken constant and the fourdimensional volume V is large, the partition 
function of gluodynamics is dominated by the contribution from the ground state, 

ZUO)=exp{-VN;^eom . (4.1) 

For Nc — »■ oo, the partition function of QCD thus approaches the expression 

Z^oM - ^ E ^M-VN^e,{d,)} , = . (4.2) 

The large V limit picks out the term in the sum for which eo{'&k) is minimal. 
This means that, in the large Nc limit, the vacuum energy density of QCD is 
determined by 

Co = min eo{^k) ■ 

k 

Since a change in 6 by 2tt is equivalent to a shift in k, the quantity Cq is periodic 
with period 2tt, while Cq is periodic only with respect to 6 ^ 6 + 27rA^. 

The distinction between the large limit of QCD and gluodynamics only 
matters for values of 6 outside the interval — tt < 6 < it. Within this region, we 
have Co = cq. For all other values, cq represents the periodic continuation from 
that interval. In our context, only the infinitesimal neighbourhood of the point 
6 = matters - we are using the external field d{x) merely as a technical device 
to analyze the properties of the Green functions for = 0. In that context, we 
do not need to distinguish between cq and cq. 

We add a remark concerning the topological susceptibility. As discussed in 
ref. P], the correlation function (0| Tuj{x) uj{0) |0) is too singular for the integral 

X{p^) = -ijd^xe'P'' {Q\Tuj{x)uj{0) |0) 

to exist, in QCD as well as in gluodynamics. The corresponding dispersive rep- 
resentation contains two subtractions, 

X(/) = X(0) +pV(0) + ^ / ,, Imx(^) . 
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Accordingly, if the susceptibility is defined as an integral over the correlation 
function, it is an ambiguous notion. We instead identify it with the mean square 
winding number per unit volume, x(0) = ('^^)/^) which does not suffer from 
such an ambiguity. In particular, if the large limit is taken at nonzero quark 
masses, we have0 

Tqcd = + 0{N-') , Tgb = 0(1) . (4.3) 

The ambiguity in the first derivative remains - it reflects the fact that, in the 
presence of a space-time-dependent vacuum angle, the Lagrangian of QCD must 
be supplemented with a contact term oc D^^OD'^O. A renormalization of the 
corresponding coupling constant is needed to absorb the quadratic divergence in 
the correlation function (0| Tuj{x) uj{0) |0). An analogous term also occurs at the 
level of the effective theory. 

5 Massless quarks 

In the present section, we briefly review some of the consequences of the U(l) 
anomaly, because these play a central role in the analysis of the low energy 
structure at large A^. For simplicity, we consider massless quarks, but leave the 
number of flavours open. We normalize the singlet axial current with 

= ^iAo7/.75g = 97^75 g , (5.1) 

where Nf is the number of flavours. The divergence of this current is given by 

= y^fU . (5.2) 

We write the 1]' matrix elements in the form 



(0| W) = tp^Fo , (0| ujW) = • 

-TO 

The relation ( |5.2|) then shows that the mass of the r]' is given by 



o 2NfT , , 

M^2, = — ^ , (5.3) 

-^0 



^The equality only holds if NcUi is large compared to the scale of the theory. In the chiral 
limit, Tqcd vanishes, so that Tqcd —/-^ t'gd- 
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and the Ward identities obeyed by the correlation functions of A° and uj lead to 



the representation 





The large Nc counting rules of section ^ imply that the correlation function of 
the axial current is a quantity of 0(i\^). For the pole term contained therein to 
be consistent with this, the decay constant Fq must be of 0[\fW^. Similarly, 
for the pole term in the correlation function formed with and lo[x) not to 

generate a contribution that diverges for ^ oo, the constant r must represent 
a quantity of 0(1). In the large limit, the pole term in (0| Ta;(a:) a;(0) |0) 
therefore disappears, as it should: That function must approach the correlation 
function of gluodynamics. In this theory, the mass gap persists when ^ oo, 
so that the expansion in powers of the momentum is an ordinary Taylor series, 
which starts with the topological susceptibility. 



In QCD, the low energy structure is more intricate, particularly when becomes 
large. In addition to the A^^ — 1 Goldstone bosons generated by the spontaneous 
breakdown of SU(Ay)flXSU(A^)i, the spectrum contains a further state that also 
becomes massless in the limit A^ — oo. At large A^, the low energy structure of 
the theory contains a new scale, M^/, which is independent of the quark masses 
and of the intrinsic scale of QCD. The magnitude of the pole terms in the above 
correlation functions, for instance, is sensitive to the relative size of the momen- 
tum p compared to this scale. 

6 Simultaneous expansion in p and 

In order to analyze the behaviour in the region where the momenta are of the 
size of M^/ - in particular, on the mass shell of the Vj - we need to treat both 
p and M^i as small quantities. This can be done in a controlled manner in the 
framework of a simultaneous expansion in powers of momenta and 1/A^, ordering 
the series with 



r = roo + 0(l/Ar,) , 



R,{j?) = 0(1) . 



(5.5) 



p = 0{V5) 



1/N, = 0{6). 



(6.1) 
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In that bookkeeping, the denominator associated with r^'-exchange represents a 
small quantity of 0{S), so that the corresponding low energy singularities are 
enhanced. In the preceding section, we considered the 1/Nc expansion at fixed p 
- as noted there, the pole term in the correlation function of a;(a;) then represents 
a subleading contribution. In the simultaneous expansion we are considering now, 
however, this term is of the same algebraic order as the one from the topological 
susceptibility of gluodynamics. Moreover, the functions Roip^) and Ri{p'^) arc 
now expanded in powers of p. Since the leading terms are of order i?o(0) = 0(1), 
Ri{0) = 0{Nc) = 0(1/6), the contributions from these functions are suppressed 
by one power of S as compared to the pole terms: For massless quarks, the 
leading term in the (5-expansion of the three correlation functions considered 
in the preceding section is obtained by simply dropping the contributions from 
Roip^) and Riip^) and is thus fully determined by the low energy constants Tqd 
and Fq. 

The result may be converted into a simple statement about the effective action. 
The correlation functions under discussion specify the part of S^ff that is bilinear 
in the external fields a|i(.T) and 9{x). At leading order in the 5-expansion, this 
part only involves the constants and Fq. The explicit expression reads 

Seff = lJdxdyf{x)f{y)A,ix-y) + lJdx{Fy^{x)a'^{x)-Tau9{xy} + ... 

where f{x) is a linear combination of the external fields d^a^^{x) and 9{x) and 
A^/(x) is the propagator of the 77': 



fix) = Fod'^alix) - ^ 0{x) , A,,{x) = j7^Jdp 



-ip-x 



Fo ^ " ' (27r)V ^ Mf^,-p^-ie ' 

The expression represents the classical action of a free scalar field -0 in the presence 
of external fields. The relevant Lagrangian is given by 

= I^D.iljDy - + e)' , D.ilj = d.ij - 2{a,) , (6.2) 



where (a^) = d^^ty^Nf is the trace of the axial external field. If we assign the 
fields the weight 

e{x) = 0(1) , a^{x) = 0{^/5) , i;{x) = 0(1) , (6.3) 

all of the contributions occurring in this Lagrangian represent terms of 0(1). 
Above, we determined the effective Lagrangian from the properties of the 

correlation functions. We could instead have investigated the general expression 
permitted by the symmetries of the theory - indeed, we will make use of that 
method to work out the higher order terms in the 5-expansion. At leading order. 
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however, that approach leads to the following problem, which arises from the 
presence of a new low energy scale and does not occur in the effective theory 
relevant for finite A^. Let us switch the external fields off. Up to total derivatives, 
the general Lorentz invariant expression quadratic in is of the form 

oo 
n=0 

A priori, it is not legitimate to dismiss terms with more than two derivatives: If 
the coefficient Cn is of order A^"', then it does contribute at the leading order of 
the 5-expansion. There is a good reason, however, why such terms do not occur. 
The point is that the propagator of the field ip is given by the inverse of the 
function ^„c„(— p^)". Unless this expression is a second order polynomial in p, 
the propagator will thus have more than one zero. With oc A^", all of these 
occur in the region where is small, of order 1/A^ - the spectrum of the theory 
would thus contain more than one particle that (i) has the quantum numbers of 
the rj' and (ii) becomes massless in the large A^ limit. Our explicit calculation of 
the pole terms in the correlation functions relies on the very plausible assumption 
that only one such particle occurs - this is why higher derivative terms do not 
occur at leading order in 6. 

7 Effective theory in the U(l) sector 

The result obtained in the preceding section shows that - as far as the two-point- 
functions of the operators A^{x) and u!{x) are concerned, and for massless quarks 
- the leading terms in the simultaneous expansion in powers of momenta and 
1/A^ can be characterized by means of a remarkably simple effective field theory. 
It involves a single dynamical variable ip{x), describing the degrees of freedom of 
the 1]'. Before generalizing this result to the other correlation functions and to 
nonzero quark masses, we add a few remarks about the structure of the effective 
Lagrangian obtained above. 

First, we note that the Lagrangian ( |6.2| ) is manifestly invariant under local 
U(A(y)jjXU(A(y)^ rotations, provided the dynamical variable is transformed ac- 
cording to 

ip' = ip -ilndetVa + iliadetVL . (7.1) 

This ensures that the sum ip + 6 is invariant, so that the same is true of the term 
proportional to Tgd- Since the trace of transforms like an Abelian gauge field, 
the term Df^ip also represents an invariant. 

Next, we observe that the choice of the dynamical variable is not unique. The 
normalization of the field ip, for instance, can be chosen such that the kinetic 
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term takes the standard form ^ d^ip d^ip (the reason for not doing so is that this 
would ruin the simphcity of the transformation law ( [7.1|) ). More generally, the 
variable if) can be replaced by a function thereof, which moreover may involve 
the derivatives of if) and the external fields - the result for the action remains the 
same. 

This implies that the effective Lagrangian is not unique. An infinitesimal 
change of variables generates a term in the Lagrangian that is proportional to 
the equation of motion. Conversely, contributions proportional to the equation 
of motion may always be removed from the Lagrangian by performing a suitable 
transformation of variables. This freedom is characteristic of effective theories, 
where the dynamical fields represent mere variables of integration in the path 
integral and are not of physical interest. The transformation ip ^ ip + n - {ip + 9), 
for instance, preserves the transformation law (|7.1| ), but takes the Lagrangian 
into 

= ^ {(1 + KfDi,^ + 2k{1 + K)DijDe + kW} - ^(1 + Kf{ij + e)\ 

with D^d = d^9 + 2(a^). For this choice of variables, the Lagrangian contains 
additional terms, proportional to D^ipD^O and D^OD^O, respectively. Note also 
that the coefficient of the one oc (^/^ + OY is then not given by the topological 
susceptibility of gluodynamics when N^. becomes large. In the following we stick 
to the form of the effective Lagrangian in eq. (|6.2| ). 

The low energy analysis relies on the perturbative expansion of the effective 
theory. At leading order, the effective action is given by the classical action, which 
collects the contributions from the tree graphs. The quantum fluctuations of the 
effective field are treated as corrections. In the standard framework, where the ry' 
does not occur among the effective degrees of freedom, the size of the quantum 
fluctuations is controlled by the external momenta and by the quark masses 
m = 0{p'^): Graphs involving i loops are suppressed by a factor of order p"^^. 
In the extended framework, however, these fluctuations involve three scales, the 
external momenta, the quark masses and the mass of the rj'. Their magnitude can 
be controlled algebraically only if M^/ is also treated as small. For the consistency 
of the effective framework, it is essential that the number of colours provides us 
with an algebraic parameter that controls the size of M^/ = 0(1/ -\/A^)- In fact, we 
will see that in the simultaneous expansion in powers of momenta, quark masses 
and 1/A^, the quantum fluctuations only start contributing at order 5^. 

The problems arising if M^/ is not treated as a small parameter can be seen 
already at tree level, where the effective action is obtained by evaluating the 
classical action at its extremum. There, the field tp{x) obeys the equation of 
motion, which is of the form □?/' + M^itp = f. At leading order in an expansion 
in powers of p at fixed M^,, this equation reduces to ip = f /M^, + . . . The kinetic 
term then only occurs as a correction. In the loop integrals, the p-expansion 
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replaces the propagators by the series (M^, —p^) ^ = + . . . It is clear that 

those integrals cannot be analyzed in this manner. 

8 Extension to U(3) 

We now discuss the extension required to analyze the full effective action and set 
Nf = 3. The dynamical variables of the effective theory must then account for all 
particles that become massless in the limit — > oo, m 0. In that limit, the 
symmetry group of the Hamiltonian is G = 11(3)^x11(3)^. We assume that this 
symmetry is spontaneously broken and that the ground state yields a nonzero 
expectation value for the quark condensate (0|gi,gfl|0). While the subgroup gen- 
erated by the vector charges remains intact the U(1)a rotations generated by 
the singlet axial charge are spontaneously broken: The operator q^qa transforms 
in a nontrivial manner under these. Accordingly, the ground state is symmetric 
only under the subgroup H = 11(3)^. The dynamical variables of the effective 
theory live on the coset space G/H = U(3), so that we may collect the effective 
fields in a matrix U{x) e U(3). The nine parameters of the coset space corre- 
spond to the nine massless pseudoscalar fields needed to describe the Goldstone 
bosons. The extension from the standard framework, where the effective field is 
an element of SU(3), to the one we are considering here shows up in the phase of 
the determinant 

det t/(x) = e^'^(^) . (8.1) 

The field tplx) describes the rj'. Under the action of the symmetry group, U{x) 
transforms with 

U'{x) = V^{x)U{x)V^{x) . (8.2) 

The transformation law ( [7. 1|) represents a special case of this formula, obtained 
by comparing the determinants of the left and right hand sides. In canonical 
coordinates, the matrix U (x) is parametrized in terms of nine pseudoscalar fields 
0O(x), ... ,08(a;): 

8 

U{x) = exp z ^ Afc0'=(x) , (8.3) 

where Ai, . . . , Ag are the Gell-Mann matrices and Aq = a/2/3. In these coordi- 
nates, the singlet field tl>{x) is represented by 'ijj{x) = v^0°(x). 

The calculation described in section ^ only concerns the part of the effective 
Lagrangian that governs the dynamics of the t]' field ipix)- The part relevant for 
the dynamics of the pseudoscalar octet may, however, be worked out in the same 
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manner. If the quark masses are set equal to zero, the eightfold way is an exact 
symmetry, so that there is no mixing between octet and singlet. The leading 
term in the 5-expansion of the correlation function of the octet components of 
the axial current, (0| Ty4^(x) ^4^(0) |0), is obtained along the same lines as in the 
case of the singlet component. The result may again be characterized by an 
effective Lagrangian that is quadratic in the corresponding dynamical variables 
- the fields (j)^{x) in this case: 

8 

^'P = lF'J2i9t.<l>'-alW(f)''-a^') . (8.4) 

k=l 

In the chiral limit, which we are considering, a mass term does not occur here: 
The pseudoscalar octet is then strictly massless. The relevant part of the effective 
Lagrangian contains a single effective coupling constant, related to the matrix 
element 

{0\A^^\n')=t6%,F . 

In fact, in the large Nc limit, the correlation functions of the singlet and octet 
currents become identical: The ratio Fq/F tends to 1. Accordingly, the sum 
C = + takes the simple form 

8 

C = lF'J2(9,<p'-a';){d>^<p'-a^'')-^{^ + ey . (8.5) 

k=0 



9 Chiral symmetry 

There is an essential difference between the Lagrangians relevant for the octet and 
singlet degrees of freedom: While is invariant under chiral transformations, 

is not. The difference arises from the fact that the part of the chiral group 
that matters for is the Abelian U(l) factor, while the one that counts for 

is a nonabelian group, 811(3)^x811(3)^. This implies that the pseudoscalar 
octet mesons necessarily interact among themselves: The Ward identities for the 
various Green functions intertwine these with one another, so that a nonvanishing 
two-point-function can occur only together with Green functions containing more 



than two currents. In contrast to the one in eq. ( |6.2|) , the Lagrangians in eqs. (B~i) 
and ( ^.51 ) cannot stand by themselves. 

The interaction terms required by chiral symmetry are well known. They 
are generated automatically, if the term 9^0'^ — is replaced by the covariant 
derivative of the effective field. For the matrix U{x) that includes the rj', the 
covariant derivative is given by 

D^U = d^U -i{v^ + a^)U + iU{v^-a^) . (9.1) 
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Under chiral rotations, this object transforms in the same manner as the field 
U{x) itself. The expansion in powers of the meson fields thus starts with 

8 
k=0 

Hence, the Lagrangian ( |8.5| ) represents the quadratic term in the chirally invariant 
expression (as usual, (A) stands for the trace of A) 

C = \F^D,U^D^^U) - ^(^ + ef . (9.2) 

The above calculation demonstrates that this Lagrangian correctly characterizes 
the two-point-functions formed with the winding number density and the singlet 
and octet components of the axial current - at leading order of the 5-expansion 
and in the chiral limit. In fact, in this limit, the Lagrangian (|9.2|) properly 
accounts for the leading contributions in the 5-expansion of all of the Green 
functions that can be formed with these operators. For a proof, we refer to 
appendix A. 

We may decompose the effective field into a part that only contains the rj' and 
a part that describes the pseudoscalar octet: 

U = e^^U, det?7 = l. (9.3) 

It is convenient to define the covariant derivative of f/ G SU(3) by 

D^U = d^U-i {v^ + a^)U + iU{v^- a^) , (9.4) 

where -0^, are the traceless parts of w^, a^. By construction, this quantity 
obeys 0^dJj) = 0, to be compared with {UW^U) = iD^^ip for U E U(3). The 
relation between the two derivatives, 

implies the identity 

which splits the above Lagrangian into an SU(3) part that exclusively involves 
the Goldstone boson octet and a U(l) part that only contains the singlet meson 
field iplx). The singlet axial field and the vacuum angle only occur in the U(l) 
part, while the vector and axial octets only appear in the SU(3) part - the singlet 
vector field does not enter at all. For massless quarks, the octet and singlet sectors 
thus separate like oil and water. 
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10 Full effective Lagrangian 



For nonzero quark masses, the properties of the theory may be analyzed in terms 
of an expansion in powers of these, so that we are then deahng with a triple 
expansion in (i) the number of derivatives, (ii) powers of quark masses and (iii) 
powers of It is convenient to generalize the ordering ( |6.1|) used for the 

massless theory by counting the quark masses as quantities of order 6: 

d^ = 0{V6), m = 0{6), 1/N, = 0{6). (10.1) 

In this bookkeeping, the vacuum angle and the effective fields U{x), iplx) are 
treated as quantities of order 1, while the external fields v^{x), a^{x), s{x), p{x) 
count as small perturbations of the same order as the derivatives and the quark 
masses, respectively: 

{U, ij, 6) = Oil) , (t;^, a^) = 0(75 ) , (s, p) = 0{5) . (10.2) 

Collecting the contributions of order 1, 5, 5^, ... , the full effective Lagrangian 
takes the form (see appendix 0): 

= £(0) + £(1) + £(2) + , _ _ (10.3) 

The individual terms only contain a finite number of effective coupling constants. 
In particular, the leading term C'^'^^ = 0(5'^) exclusively involves the coupling 
constants F, B and r. Replacing the external fields s{x) and p{x) by 

x{x)=2B{s{x) + ip{x)} , (10.4) 

the explicit expression reads 

£(0) ^ l^F^^D^U^D^U) + \F\U^x + X^U) - | r(^ + Of . (10.5) 

Note that C'^'^^ contains vertices of the type F'^dcjydcjxjf^''^ ^ which represent in- 
teractions among the pseudoscalar mesons. The corresponding tree graph contri- 
bution to the scattering amplitude describes a collision with altogether n particles 
in the initial and final states. In the large Nc limit, the contribution disappears, in 
proportion to oc Nc~^^'^, in accordance with the general counting of powers 
for the scattering amplitude given in appendix 0. When is sent to infinity, all 
of the bound states become free particles of zero width. The larger the number of 
particles participating in the reaction, the smaller the interaction among these. 

In the large limit, diagrams with a single quark loop dominate, so that only 
contributions with a single trace over flavour matrices occur (OZI rule). This 
property entails that the two-point-function of the singlet current approaches the 
one of the octet components, in particular Fq/F — > 1. In general, however, the 
matrix elements of the singlet and octet components are different, even if Nc 
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is taken large. In particular, the OZI rule does not imply that the four-point- 
function of the singlet axial current approaches the one of the octet components 
- in fact, it does not: The latter contains a four-fold pole, whose residue is 
proportional to the scattering amplitude. The term represents a contribution of 
order A^, which at low energies is fully determined by the pion decay constant. 
An analogous contribution to the four-point-function of the singlet axial current 
does not occur, because the scattering amplitude T^'r/' rj'r]' does not pick up a 
term proportional to Crossing symmetry implies that, for massless quarks, 

the coefficient is proportional to s + t + u = 4M^, = 0(1/A^). 

The term C^^^ = 0{S) contains the contributions of 0{Ncp'^), 0{p^) and 
0{1/Nc). The first may be copied from the SU(3) Lagrangian, simply dropping 
those coupling constants that violate the OZI-rule. Concerning the remaining 
contributions, the terms that are of the same structure as those in C^^'^ may be 
absorbed in the coupling constants F,B,t (note that the coefficient of the term 
proportional to [ijj + 6)"^ then differs from the topological susceptibility of gluo- 
dynamics - we need to distinguish the coupling constant r from Tgd)- With a 
suitable choice of the dynamical variables, the explicit expression for C^^^ can be 
brought to the form 

£(^) = L2{D^U^DMD^'U^D''U) + (2L2 + L^){D^U^ D^'U D^U^ D^U) 
+L,{D^U^D>^U{U^X + X^U)) + LsiU^xU^X + X^Ux^U) 
-iLg{R^,D^UD''W + L^^D^U^D'^U) + L^^^R^.U L^^'U^) (10.6) 
+ ^F^AiD^iPD^^ - + e){U^x - X^U) 

+ l^HoD^eD^e + H^iR^^R^'' + L^^Ln + H^ix^x) ■ 

The covariant derivatives and the field strength tensors are defined by 

Df,U = d^U - i {Vf, + a^)U + iU{Vf,- a^) , 

D^^l, = d^^ - 2(a^) , D^e = d^e + 2{a^) , (10.7) 

Ri^u d^Tj^ dj^T^ ^[^/i; 7 L^i^ d^li/ dylpi_ i\j'fj,i lu\ > 

where = + and = — a^. The somewhat queer numbering of the 
coupling constants arises because we retain the notation introduced in ref. |^ to 
denote the couplings at first nonleading order. The SU(3) Lagrangian given in 
that reference contains three independent terms with four derivatives, involving 
the coupling constants Li, L2 and L3, respectively. For 00, all of these 

are of order A^, but the OZI rule suppresses the combination 2Li — L2. The 
above expression for those terms that involve four derivatives is obtained from 
the Lagrangian of ref. by replacing the constant Li by \L2. The counting 
rules imply that the coupling constants L2, L^, L^, Lq, Lg, Liq represent quantities 
of O(A^), while Ai, A2 are of 0(1/A^). Concerning the contact terms, Hq is of 
order 1, while Hi and H2 are of O(A^). 
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The number of independent coupling constants entering the low energy rep- 
resentation of the effective action to first nonleading order is roughly the same 
as for the SU(3) Lagrangian: While the four parameters Li — |L2, L4, Le, are 
relegated to next-to-next-to leading order, three new coupling constants appear: 
r, Ai,A2. The framework allows us to also evaluate the correlation functions 
formed with the operators ^4° and tu, while the Lagrangian of ref. accounts 
for these only at leading order. A further virtue of a Lagrangian that explicitly 
includes the degrees of freedom of the rj' is that it equips the contributions from 
r^'-exchange with the proper denominator and, for instance, distinguishes the fac- 
tor 1/M2, from 1/(M2, - M^). In view of M^/M^, = 0.33, this distinction is 
numerically quite significant. On the other hand, the above framework retains 
only the next-to-leading terms of the expansion. In sections |1^ and 



we 



will consider those contributions of next-to-next-to leading order that are needed 
to fully match the SU(3) and U(3) theories on their common domain of validity. 



11 Potentials 

At this point, we establish contact with the notation introduced in ref. 0], where 
the vertices of the effective Lagrangian are ordered in the same fashion as for the 
case of SU(3), that is by simply counting powers of momenta and quark masses. 
We refer to that ordering as the p-expansion. Ordering the various couplings in 
this manner, the Lagrangian takes the form 

Ceff = -Vo + Vi (D.U^D^'U) + V2 (U^x) + ix^U) (11.1) 

The leading term Vq contains those vertices that do not involve derivatives or 
quark mass factors. Chiral symmetry implies that the matrix U cannot appear 
without derivatives, but it does not protect the combination ip + 6 oi the singlet 
field and the vacuum angle. There are vertices without derivatives containing an 
arbitrary power of this combination. Their collection Vq represents a function of 
the variable ip+9. The same applies for the coefficients of those contributions that 
do involve quark mass factors or derivatives. The effective coupling constants are 
the coefficients Ki,fc occurring in their expansion in powers oi ip + 9, 

00 

k=0 

The functions Vn{ip + 9) may be viewed as potentials that describe the dynamics 
in the U(l) sector, with Vq as the most important one. The counting rules 
established in appendix ^ imply that - for a suitable choice of the dynamical 
variables - we have 



19 



Indeed, these rules were written down already in ref. [^. The detailed analysis 
described in the present paper merely fills a gap: The argumentation given in that 
reference is incomplete, because it does not take the ambiguities into account that 
arise from the freedom in the choice of the dynamical variables. As demonstrated 
in the appendix, these ambiguities are in one-to-one correspondence with those 
occurring in the off-shell extrapolation of the matrix elements. What we have 
shown now is that there exists a class of such extrapolations, for which the above 
counting rules for the effective coupling constants strictly follow from those for 
the correlation functions. We will make use of this result when analyzing the 
Kaplan-Manohar transformation. 

At leading order, the freedom in the choice of the dynamical variables is hid- 
den in the manner in which the field U{x) is parametrized in terms of the variables 
. . . , - that is in the choice of the coordinates on the group U(3). As 
we are expressing the Lagrangian in terms of U{x), it is irrelevant how these coor- 
dinates are chosen. At first nonleading order, however, the ambiguities inherent 
in the choice of the dynamical variables do show up and the Lagrangian takes a 
unique form only if that freedom is fixed. As discussed in appendix |^, this can 
be done by eliminating invariants that vanish if the field U{x) obeys the equation 
of motion associated with C^^^ . The representation for C^^^ in eq. ( |10.6| ) fully ex- 
ploits this freedom. The eleven couplings occurring there are both complete and 
independent: With a suitable change of variables, all of the invariants that can 
be built up to and including 0{6) reduce to a linear combination of those listed 
and, conversely, the choice made fixes the dynamical variables - up to transfor- 
mations of 0(5^). In the notation introduced above, the choice of variables made 
in eq. ( |lU.(i| ) corresponds to V4 = 0, as far as the terms of 0{6) are concerned. 



One readily checks that a suitable change of variables of the type U — > e'^^U, 
where / only depends on ip + 6, removes V4 to all orders and thereby fixes the 
dynamical variables up to transformations of 0{p^). 

At order p^, the natural parity part of the effective Lagrangian contains al- 



together 52 potentials. For an explicit representation, we refer to [T^]. The 
potentials relevant for the unnatural parity part are given explicitly in section |2T. 
In connection with the Kaplan-Manohar transformation, the symmetry breaking 
terms of order are of special interest. These are of the same structure as in 
the case of the SU(3) Lagrangian: 

C,2 = WiiU^xU^x) + W*{x^Ux^U) + W2{U^x)' (11.2) 
+W^{x^U)^ + Ws{U^x){x^U) + W,{xh) . 

The only difference is that the coupling constants Lq,Lj, Lg and H2 are replaced 
by the potentials Wi, . . . , W4, which depend on the variable ^ip + 9. 

In the expansion, most of the potentials occurring in C2 start contributing 
only at next-to-next-to leading order: If only terms of order 6 are retained, the 
effective Lagrangian reduces to the explicit expression in eq. ( |10.6| ). 
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12 Dependence on the vacuum angle 6 



Since the dependence of the various quantities on the vacuum angle is rather 
pecuhar, particularly in the large limit, we now briefly discuss this issue. 
Throughout this section we restrict ourselves to a constant value of the external 
field 6{x) = 9, turn the pseudoscalar field p{x) off and identify the scalar one 
with the quark mass matrix, s{x) = m. 

The transformation law (|2.3| ) shows that a suitable global U(l) transformation, 
for instance Vr = e^^l, = 1, removes the vacuum angle. The external vector 
and axial fields stay put, but the quark mass matrix undergoes a change: m — > 
mg. The invariance of the effective action thus implies 

S{v, a, 6, m} = S{v, a, 0, mg} , mg = e^^m , (12.1) 

the familiar statement that the vacuum angle only enters in combination with 
the quark mass matrix: Only the phase arg det mg = arg det m + 9 matters. 

In the effective Lagrangian, however, the vacuum angle does not enter in 
this combination. Also, according to eq. ( ^.2|) , the effective field undergoes a 
change under the above transformation: U — > e^^U. We need to be careful when 
performing finite transformations on the effective field. In the construction of 
the effective Lagrangian, only the series in powers of both the external fields and 
the dynamical variables (pix) were considered. If we now take a vacuum angle of 
finite size, we are summing the expansion in powers of 6 to all orders. This may 
lead to ambiguities. The relation detU = e*^, for instance, can be solved for ip 
only up to a multiple of 2it. This means that the decomposition in eq. ( |9.^j| ) is 
unique only in the vicinity of U = 1. 

Since the effective theory is evaluated perturbatively, its properties are gov- 
erned by the tree graphs, that is by those of the corresponding classical field 
theory. In this framework, ambiguities of the type just mentioned are avoided 
if is replaced by the gauge invariant field ip = ip + 6, collecting the remaining 
dynamical variables in the matrix U : 

U = e^^U, ip = ilj + e, detU = e-'^. (12.2) 
It is convenient to define the covariant derivative of U by 

D^U = d^U-t{v^ + a^)U + iU{v^-d^) , (12.3) 

SO that (JJ'^D^U) = 0. In this notation, the derivative of U reads 

D,U = e^^ {D,U + iid^tp - D^e)U] , (12.4) 
and the leading term in the effective Lagrangian takes the form 
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The equation of motion belonging to this Lagrangian unambiguously determines 
the classical solution ipd, Ud in terms of the external fields. If these are subject to 
a gauge transformation, ipd remains invariant, while Ud goes into V^f/czVcf, even 
if the angles of the transformation are taken large. This also holds if the higher 
order contributions to the effective Lagrangian are accounted for. In particular, 
under the global U(l) rotation considered above, which maps the vacuum angle 
into zero, the classical solution does transform with Ud — > e^^Ud- 

At leading order of the perturbative expansion, the effective action is given by 
the value of the classical action at the extremum. Hence it does obey the relation 
( |12.1 ), also if 9 is taken large. Note that, once the dynamical variables of the 



effective theory are identified with ip and U , the effective Lagrangian contains 
the external field 6{x) exclusively through the derivatives thereof - a constant 
vacuum angle then exclusively manifests itself via the constraint det f/ = e~*^. 
This also implies that the effective theory automatically generates an effective 
action that is periodic in the vacuum angle: The only place where a constant 
shift in that angle shows up is through the factor e~*^, which remains the same 
if 9 is replaced by + 2tx. 

The dependence on 9 disappears in the large limit, both in gluodynam- 
ics and in QCD. The vacuum energy density, for instance, represents a term of 
0{N^), but a dependence on 9 only shows up at 0(1), through the susceptibility 
term. The origin of this suppression can immediately be seen in the structure 
of the QCD Lagrangian: The properties of the theory are governed by the term 
g~'^G^vG^'^ , while the vacuum angle enters through 9G^vG^^. In the large Nc 
limit, the term containing the vacuum angle is smaller than the one that deter- 
mines the dynamics by the factor g'^9 = 0{1/N^. 

Gauge invariance implies that, in the effective Lagrangian, the vacuum angle 
only appears in the combination ip = ip + 9, together with the singlet field. In 
particular, the term proportional to 9"^ in the vacuum energy density is converted 
into one proportional to ip"^ and thus equips the r]' with a mass. That mass 
disappears in the large limit, because the 6'-dependence of the vacuum energy 
of gluodynamics is suppressed. 

Small quark masses also suppress the dependence on 6*. If the determinant 
of the quark mass matrix vanishes, the Green functions of the vector and axial 
currents and hence also the scattering matrix elements even become entirely inde- 
pendent of 9. Goncerning the 6'-dependence, the quark masses - the parameters 
that break chiral symmetry - thus play a similar role as the parameter l/N^ that 
measures the breaking of the OZI rule. 

As an immediate corollary of the fact that, in the large limit, the depen- 
dence on the vacuum angle is suppressed, the potentials reduce to constants in 
that limit. More generally, if the expansion in 1/A^ is cut off at a finite order, 
the potentials are replaced by polynomials in ip. 
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13 Kaplan- Manohar transformation at large 



As pointed out by Kaplan and Manohar the standard SU(3) Lagrangian is 
invariant under the transformation in eq. ( |1 . 1| ) , provided the effective couphng 
constants Lq, Lj and Lg are subject to a corresponding change. In that framework, 
the vacuum angle is set to zero. The transformation can, however, be generalized 
to allow for nonzero values of 9{x). To ensure that the modified mass matrix 
has the same transformation properties as the original one also with respect to 
chiral U(l)-rotations, it suffices to equip the transformation law for the scalar 
and pseudoscalar external fields with a factor of e~'^. With m{x) = s{x) +ip{x), 
the Kaplan-Manohar transformation then takes the form 

m — m + A e~*^ m^~"^ det . (13.1) 

In the present context, it is convenient to order the vertices with the p- 
expansion, using the representation of the effective Lagrangian in eq. At 



order p"^, the scalar and pseudoscalar external fields only enter trough the term 
y2{U^x) ciiid its complex conjugate, with x = '^Bm. Under the above operation, 
this term picks up a contribution proportional to det ■ The expression 

may be simplified by applying the identity 



2 detC 

to the matrix C = x'^U. In view of det U = e^^, the transformation then yields 

V2{U^X) - V2{U^x) + ^V2e-'^^^'^ {{x^Uf - {xWU)} . 

The modification may be absorbed in a suitable change of the potentials that de- 
scribe the mass terms contained in the effective Lagrangian of order p^, which are 
given in eq. ( |11.2D . Indeed, the effective Lagrangian does remain invariant under 
the transformation ( |13.1| ), provided the effective coupling constants contained in 
the potentials Wi and W2 are modified according to 

Wi^Wi + ^ e*('^+^) , (13.2) 
AB 

Aid 

all other coupling constants staying put. This demonstrates that the general 
U(3) Lagrangian exhibits the same kind of reparametrization invariance as the 
one built on SU(3). 

The modification of the potentials Wi, W2, however, is in conflict with the 
behaviour of the theory at large As discussed in detail in the preceding 
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section, the effective Lagrangian represents a polynomial in the variable ip + 6, 
at any order of the expansion - an immediate consequence of the fact that 
the dependence of the matrix elements on the vacuum angle is suppressed in the 
large Nc limit. In view of the factor e^^'^~^^\ the changes required in Wi and W2 
violate this condition, already at leading order. The disaster evidently originates 
in the phase factor e~*^ occurring in the Kaplan-Manohar transformation ( |13.1| ): 
This factor is needed for the modified quark mass matrix to transform in the 
same manner as the physical one, but it introduces a dependence on 6 that is in 
conflict with the large Nc properties of the theory. 

It does not come as a surprise that the Kaplan-Manohar transformation breaks 
the Okubo-Iizuka-Zweig rule: The transformation mixes different quark flavours. 
It is to be expected that the parameter A can at most be of order This also 

follows from the fact that the coupling constant Lq picks up a term proportional 
to XF"^. Since the counting rules imply that Lg represents a term of 0(1), while F"^ 
is of order A^, the transformation can be consistent with the large properties of 
the theory only if A represents a term of order The above analysis reaches 

much further: It shows that the parameter A vanishes to all orders of the 
expansion. 



14 Renormalization of the operators and lj 



As is well-known, the dimension of the singlet axial current is anomalous p9 
The operator ^4° = q^Xo'y^'j5q must be renormalized for the correlation functions 
formed with this current to remain finite when the cutoff is removed. In particu- 
lar, the one particle matrix elements (0| A^^\P) = ip^Fp require renormalization: 
The singlet decay constants F°o, represent quantities like the quark masses 

or the quark condensate - they must be renormalized. 

In the MS scheme, the renormalized quantities depend on a running scale. 
To distinguish this scale from the one used to renormalize the loop graphs of 
the effective theory, we denote it by /Iqcd- There is a difference in that the 
scale dependence of the quark masses shows up already at leading order in the 
1/A^ expansion, while the anomalous dimension of the singlet axial current only 
manifests itself at next-to-leading order, because the triangle graph responsible 
for the phenomenon contains an extra quark loop. 

The representation of the renormalization group only interrelates operators 
with the same dimension and the same Lorentz quantum numbers. Moreover, 
the matrix elements of the representation exclusively involve those coupling con- 
stants that are dimensionless in the dimension of physical interest - the QCD 
coupling constant g and the vacuum angle 9 in the present case. In particular, 
the representation is independent of the quark masses. If these are turned off, 
the charges of the flavour group SU(3)ijxSU(3)i are conserved. The correspond- 
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ing representation on the set of field operators commutes with the one of the 
renormahzation group. 

Since the axial singlet current does not have any partners with the same 
Lorentz and flavour quantum numbers and the same dimension, it transforms 
as a singlet under the renormahzation group. The operators g^g^ form an irre- 
ducible representation of the flavour group. As the theory does not contain any 
other operators with these quantum numbers and with the same dimension, the 
renormahzation is multiplicative also in this case, 

= Z^Al , (g^g^)-" = Z,,g^g^ . (14.1) 

The renormahzation of the singlet decay constants, for instance, reads 

The renormahzation factors depend on the running scale of QCD: 

,^.-;!^-..Z., -,,._5^M^(A)Vo(,V ,14.2) 

d/iQCD Nc V47r/ 

^'^^^^ = ' = + ) ■ ^^'-'^ 

The renormahzation of the operator uo is more complicated. As this field 
represents the variable conjugate to 6, the issue is related to the dependence of 
the effective action on the vacuum angle. This dependence is of crucial importance 
for the analysis of the theory at large A^. For this reason, we now discuss it in 
detail. 

The external fields may be viewed as space-time dependent coupling con- 
stants. The renormahzation of the operators amounts to a renormahzation of 
these "constants" . In particular, if the vacuum angle is turned off, the effective 
action remains the same if we replace the bare operators g^g^, g^g| and A° by 
the renormalized ones and at the same time replace the bare external fields by 
the quantities s^*^"" = Z~^s, p^*^"" = Z~^p, {a^y^"^ = Z^^(a^), while the vector 
field and the octet components of the axial field remain put - the corresponding 
operators do not get renormalized. 

To see what happens if the vacuum angle 6 is turned on, we exploit the fact 
that the natural parity part of the effective action is invariant under the local 
U(3)flXU(3)i transformations specified in eq. ( p.3|) . We denote the octet part of 
the axial field by = — and replace s, p by the combination 

= e^^(s + ip) , (14.4) 
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which is invariant under the transformations generated by the singlet axial charge. 
Invariance under this subgroup then implies that (a^) and 9 can enter the effective 
action only through the gauge invariant combination D^O — d^j^O + 2(a^), 

Seff = Seff{v^, a^, D^e, mj (14.5) 

This shows that the functional Seff is known for arbitrary 9{x) if it is known for 
9{x) = 0: The general expression is obtained from the one relevant for 9{x) = 
by replacing (a^), s and p with \D^9, \ij^e + ^^"^ f (^J ~ "^^e)' respectively. 
In particular, the quantities 0^9 and me are renormalized according to 

{D,9)- = Z-'D,9 , m;- = Z'/m, . (14.6) 

15 Renormalization of the effective action 

For the effective action to remain finite when the cutoff is removed, the Lagrangian 
must inchidc all terms of mass dimension less than or equal to four that are 
consistent with the symmetries of the theory. This also applies to contact terms, 
such as hQD^9D^'9, hi{R^^R>''' + L^^L^'^) or h2{m)m). The constant hi, for 
instance, is needed to renormalize the QCD contributions to the electric charge. 
It generates contact terms in the correlation functions of the vector and axial 
currents. Together with the renormalization of the bare coupling constant the 
operation ensures that all of the renormalized correlation functions formed with 
the quark currents and with uj approach finite limits, 

SeffW, a„ D,9, m„ g, h, /xo} = 5e//{^^, a„ {D,9Y'\ m^", //qcd}- (15.1) 

We have denoted the coupling constants associated with the contact terms col- 
lectively by h and use the MS scheme for the renormalized quantities. These are 
independent of the cutoff /xq, but do depend on the running scale /Xqcd- In the 
effective action, the running scale drops out - the left hand side of the above 
relation is scale independent. 

A priori, since 9 is dimensionless, the representation of the renormalization 
group might depend on this coupling constant, too. As discussed above, how- 
ever, 9 is an inessential parameter in the sense that the manner in which the 
effective action depends on it is fully determined by the symmetries of the the- 
ory. In particular, the factors and Z^^ arc the same as for = 0. In fact, the 
renormalization law for rriQ shows that the vacuum angle is a scale independent 
quantity - in agreement with the fact that the effective action is periodic in this 
variable. 

The renormalization law for D^9, on the other hand, implies that the renor- 
malization of the singlet axial field is rather complicated: The effective action 
approaches a finite limit only if, in addition to a multiplicative renormalization 
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with the factor ^, the field (a^) is simultaneously also subject to a U(l) gauge 
transformation by the angle ^{Z~^ — l)9{x), 

{a,r- = Z-\a,) + i - 1) d,e . (15.2) 

The relation concisely specifies the renormalization of the correlation functions 
involving the operator u. In contrast to and g^g^, for which the renormaliza- 
tion is multiplicative, uj is subject to an inhomogeneous renormalization, despite 
the fact that the variable conjugate to it, the vacuum angle, is renormalization 
group invariant. 

The origin of the complication is readily understood: It reflects the conserva- 
tion law of the singlet axial current 

9Mj = (7° + VGuj , (7° = sjlqil^mq (15.3) 

(for simplicity, we identify s with the quark mass matrix m and set p = 0). While 
d^A^^ picks up a factor of Z_^, the term is invariant, because the qiiark mass 
matrix transforms contragrediently to the operator qi'j^q. The matrix element 
{0\uj\ri'), for instance, may be represented as a difference of two contributions. 
The first is proportional to and thus scales with Z4, while the second is 

given by the matrix element (0| a \r]'), which does not get renormalized. 

In contrast to the operator uj itself, the zero momentum projection thereof, the 
winding number v = Jdxu, is renormalization group invariant. The correlation 
functions of this quantity may be obtained from the effective action by considering 
a constant vacuum angle and taking derivatives with respect to it. The term 
proportional to d/j^O occurring in the renormalization of (a^) is then absent. Since 
the vacuum angle is invariant under renormalization, the same holds for the 
correlation functions of p. 

To illustrate the statement, consider the topological susceptibility of QCD 

Tqcd = -iJdx{0\Tu;{x)u;{0) |0) . 

If s is identified with the quark mass matrix, the vacuum angle is taken con- 
stant and all other external fields are switched off, the effective action reduces to 
Seff = — J dxe{m,9), where e{m,9) is the vacuum energy density of QCD. The 
topological susceptibility is the second derivative of e(m, 9) with respect to 9. It 
can be represented in terms of the correlation function {0\Ta^{x) a^{y) |0) and 
the quark condensate, by invoking the Ward identities 

9^(0|TAj(x)cu(|/)|0) = {0\Ta\x)Luiy)\0) + V6{0\TLoix)Luiy)\0) 

d''{Q\TAl{x)a'{y)\0) = {0\Ta%x) a%y) \0) + V6 {0\Tujix) a\y) \0) 

—i |5(a; — y){0\qmq |0) . 
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Since the left hand sides represent total derivatives, they drop out when taking 
the integral over all of space. This leads to the representation 



r,co = ^ JdxiO] Tcuix) a%0) |0) = jdx{Q\ Ta%x) a%0) |0) - ^ (0|gmg |0). 

The relation confirms the statement that the susceptibility is renormalization 
group invariant: Neither the correlation function (0| Ta^{x) cr*^(0) |0) nor the term 
involving the quark condensate, {0\qmq |0), depend on the running scale of QCD. 

16 Dependence of the effective coupling con- 
stants on the running scale of QCD 

We now translate these properties of the effective action of QCD into the lan- 
guage of the effective theory. The fact that the operators q^q'l, A° and u must 
be renormalized implies that some of the effective coupling constants depend on 
the running scale of QCD. Apart from contact terms, the renormalization is fully 
determined by the factors and Z^^ - the coupling constant g is hidden in the 
effective couplings. The contact terms ho, hi,... of the QCD Lagrangian are 
absorbed in the coupling constants Ho, Hi, . . . of the effective theory. The renor- 
malization of these couplings thus involves the renormalization factors relevant 
for ho,hi, . . . In the following, we disregard the contact terms altogether - they 
lead a life of their own. Note that the present section concerns the scale depen- 
dence of the effective coupling constants that arises from the renormalization of 
the QCD Lagrangian. The one generated by the logarithmic divergences that 
occur within the effective theory is an entirely different issue (see section p!7|). 

If the external fields (a^) and 9 are switched off, the renormalization exclu- 
sively concerns s and p. For the generating functional of the effective theory 
relevant at fixed to become renormalization group invariant, the low energy 
constant B must be renormalized with = Z^^B. The scale dependence of 
this constant cancels the one of the fields s{x) and p{x), so that the quantity 
X = 2B{s + ip) is independent of the QCD scale. If the effective SU(3) La- 
grangian is written in terms of this variable [Q, all of the coupling constants 
occurring therein, the pion decay constant in particular, are invariant under the 
renormalization group. 

Let us now consider the extension to U(3) and include the external fields (a^) 
and 6, so that the natural parity part of the effective action becomes invariant 
under local U(3)flXU(3)^. The renormalization properties discussed in the pre- 
ceding section do not rely on the large limit. It is therefore appropriate to 
return to the form of the effective Lagrangian in eq. ( |11.1| ), which does not invoke 
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the expansion. Inserting the decomposition ( 



123), we obtain (V^ = 0): 



+ilv^ + v,){d^i^ - D^ef + v,D^eD>^e + o(p^) . (le.i) 

The covariant derivative -D^t/ involves the octet components of the vector and 
axial fields, as well as the vacuum angle. Neither one of these quantities under- 
goes renormalization, but the singlet axial field in D^6 = 8^9 + 2(a^) does: In 
addition to a multiphcative renormalization, it picks up a U(l) gauge transfor- 
mation proportional to 9^^, so that [D^Oy^'^ = Z^^Dfj^O. For the operation to 
map the effective Lagrangian onto itself, the variable ip must be renormalized in 
the same manner: 

= Z-'ij , (16.2) 

while U remains put. The Lagrangian then remains invariant, provided the ef- 
fective coupling constants are renormalized in such a manner that the potentials 
transform with 



V.ixY'- = V,{Z,x) 

V2{xY''' = V2iZ^x)e~^^^-^-'^'' (16.3) 
VsixY'^ = ZlV,{Z^x) + \{Zl-l)Vr{Z^x) 
V^ixY''' = ZlV,{Z^x) . 

We conclude that the effective Lagrangian is invariant under the renormaliza- 
tion of the external fields specified in section |1^, provided the dynamical variable 
U is subject to renormalization: While the unimodular part remains invariant, 
the phase det U = e^"^ transforms with 

Expressed in terms of U, the renormalization thus amounts to 

Ijren ^ ^^^^ Uf^ U , Z = - 1 . (16.5) 

We repeat that we are disregarding contact terms. The above effective La- 
grangian does contain such a term: Hq = 12V5(0). The renormalization of Hq 
involves the one of the QCD counter term ho and is not covered by the above 
renormalization prescription for the potential V5. Except for Hq, the relations 
( |16.3|) do specify the renormalizations of all of the effective coupling constants 



occurring to order p^, in terms of the factors and Z^^ that characterize the 
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anomalous dimensions of the operators A° and g^g^. The couphngs collected in 
the potential Vq, for instance, are renormalized according to 

For the constant r = 2 V02, this yields 

Note that the factor Za differs from 1 only by a term of order As the 

triangle graph responsible for Z^ does not occur in gluodynamics, Tgd does not 
get renormalized. 

The transformation property of Vi = jF'^+0{x'^) shows that F is independent 
of the QCD scale - in agreement with the fact that F represents the value of a 
physical quantity in the chiral limit. The constant B transforms contragrediently 
to the quark mass matrix, B"^^^ = Z^^B, so that the lowest order mass formula 
= {niu + md)B + . . . does yield a scale independent pion mass. For the 
coupling constants Ai and A2 the above relations yield 

l + Ar = Zi(l + Ai) , l + Ar = Z,(l + A2) . (16.7) 

The result shows that the renormalization of the effective coupling constants is in 
general not multiplicative. This reflects the fact that the individual terms of the 
series C^^^ + C^^^ + . . . are not invariant under the renormalization group. The 
leading term for instance contains the contribution j^F'^Df^ipD^ip, which picks 
up a factor of Z'"^. A term proportional to D^ipD^'ip also occurs in C^^\ The 
two contributions add up to + KijO^tpD^ilj - the renormalization of Ai 

indeed ensures that the sum is renormalization group invariant. Likewise, the 
renormalization of the potential V2 generates a term proportional to Z^^ — 1, 
which is absorbed in the renormalization of A2. The renormalization group thus 
intertwines terms in the effective Lagrangian that involve the same number of 
derivatives and quark mass factors, but carry a different power of A^. This should 
barely come as a surprise - the renormalization factor Z^ itself represents an effect 
of the type 1 + 0(1/A^). The renormalization of the nonleading couplings is not 
multiplicative, because it must cure the deficiencies of the leading terms, so that 
the results obtained on the basis of the first two terms in the 1/Nc expansion 
of the effective Lagrangian do become renormalization group invariant to first 
nonleading order. 

The remaining coupling constants of C^^^ are renormalization group invariant, 
because the renormalization of C'^'^^ does not generate terms of order p^. They do, 
however, give rise to specific nonmultiplicative renormalizations of the couplings 
occurrmg m 

£{2) 

(see section |l^). 

Strictly speaking, the preceding discussion only shows that the renormaliza- 
tion of the effective coupling constants which we have just given is sufficient for 
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the effective action to be independent of the QCD scale. One may also show that 
this renormalization is necessary - it suffices to calculate a few observables within 
the effective theory. The scale independence of the result indeed implies that the 
effective coupling constants must be renormalized in the above fashion |jl9|. 



17 Higher orders and loops 

There is a significant difference between the U(3) framework considered in the 
present paper and the standard one, where the degrees of freedom of the meson 
field are restricted to those of SU(3). In that case, the loop graphs of the effective 
theory are relevant already at first nonleading order. Now, they only matter if we 
wish to extend the calculation beyond this order. As far as powers of momenta 
are concerned, the series is of the same type in the two cases: Momenta count 
like p ~ The loop graphs, however, are inversely proportional to powers of 
~ While in the standard chiral perturbation series, graphs containing i 

loops generate contributions of order p"^^, they now only manifest themselves at 
order 5^^. In particular, the one loop graphs yield contributions of next-to-next-to 
leading order. 

Numerically, the one loop graphs are of about the same size as those occurring 
in the chiral perturbation series of SU(3) - to count F^^ as a term of order a/A^ 
does not change the numerical value of this coupling constant. The one loop 
graphs are relevant also if we wish to establish the relation between the coupling 
constants occurring in the two versions of the effective theory: The SU(3) coupling 
constants depend on the running scale used to renormalize these - we cannot 
make a significant comparison if we ignore the loops. The evaluation of the 
one loop graphs, however, poses a problem: Since these violate the OZI rule, 
the divergences can be absorbed in a renormalization of the effective coupling 
constants only if we include terms of 0{p'^) that are sub leading in 1/Nc. 

A complete analysis of the contributions of next-to-next-to leading order is 
beyond our scope. The corresponding part of the effective Lagrangian contains 
a plethora of terms of order NcP^, p^, p^/N^ and respectively. The first 

category contains vertices of the same structure as those occurring in the general 
SU(3) Lagrangian of order p^, which are listed in ref. For a complete set 

of terms of order p^, we refer to As the one loop graphs of C^^^ represent 



contributions that involve at most four powers of momentum, their renormaliza- 
tion only requires counter terms that are at most of order p^. In other words, the 
coupling constants of the type NcP^ are independent of the running scale of the 
effective theory. Also, since the QCD renormalization group preserves the Lorentz 
structure of the various vertices, it can intertwine the couplings of C^^^ and C^^^ 
only with those terms in C^'^^ that are at most of order p^. In the following, we 
restrict ourselves to a small subset of the couplings occurring at next-to-next-to 
leading order, the minimal set needed to satisfy the following requirements: 
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(a) All of the couplings relevant for those quantities that can be calculated 
within the framework of ref. |Q are included. This allows us to match the two 
versions of the effective theory at one loop level and to express all of the SU(3) 
coupling constants in terms of those occurring in the U(3) Lagrangian. 

(b) All contributions needed to absorb the divergences occurring in the one 
loop graphs for the masses and decay constants of the pseudoscalar nonet are 
included, so that we can unambiguously account for those contributions of next- 
to-next-to leading order that are enhanced by a chiral logarithm. The terms 
required by (a) suffice to satisfy this condition, except for one extra coupling. We 
retain the numbering introduced in ref. Il^ and denote this term by Lig. 



(c) The conditions (a) and (b) ensure that the result for the decay constants 
is invariant under the renormalization group of QCD (see section |T6|). For this 
to be the case, it is essential that the coupling constant Lig is included: The 
renormalization group intertwines L5 with this term. Likewise it intertwines Lg 
with a further extra coupling constant, L25, which needs to be included for the 
masses of the pseudoscalars to also become renormalization invariant. 

The effective Lagrangian then contains the following couplings: 

£(2) = {U-\L2){D^U^D^Uf + L,{D^U^D^U){U^X + X^U) 

+L,{U^X + X^Uf + L^{U^x - X^Uf (17.1) 
+Lig iD^ilj {D>^U\ - D^Ux^) + L25 ^(^ + &) {U^xU\ - X^Ux^U) . 

As all of these terms represent contributions of O(p^), the coefficients approach 
a finite limit when Nc —>■ 00. The contributions involving Li — |L2, L4, Lg and L7 
correspond to those pieces of the SU(3) Lagrangian that violate the OZI-rule. 

The divergences generated by the one loop graphs of C^^^ have been worked 
out in ref. fl^. Expressed in terms of the factor 



{-^-i(ln4. + r(l) + l)} (17.2) 



167r2 [d-4: 

the renormalization of the effective coupling constants required to absorb these 
divergences takes the form: 

5 = 5^|l + , Ln = Ll + r„A , Hn = K + A„A . (17.3) 

The constants F, r, Ai and A2 do not get renormalized. The main difference 
to the case of SU(3) is that loops involving the propagation of an 77' require a 
renormalization of the low energy constant B. The same graphs also generate a 
change in the values of the coefficients Fe, Fg and A2: 

Ti = F2 = A, F3 = 0, F4 = F5 = |, Fg = TTT, F7 = 0, Fg = -j^. 



32' ^ 16 ' •^ "5-4 g5^o 8' '^ 16' ' 5 

i? Tio = — |, Fig = — i, F25 = 0, Ao = 0, Ai = — |, A2 = I . 
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These relations concern the renormahzation of the effective couphng constants 
required to absorb the divergences occurring within the effective theory - the 
dependence on the renormahzation scale of QCD is a different matter. For the 
coupling constants of C^^^ and C^^\ we have discussed the issue in detail in section 
That analysis is readily extended to the terms contained in C'^'^\ We only 
give the result: (a) The coupling constants Li,L4,Lq and L7 are independent of 
the QCD scale, also in the framework of U(3). (b) The constants Lig and L25 
must be renormalized according to 

2L5 + 3Lll^ = Z^{2L, + 3L,s) , 2Ls - SL^^ = W^Ls - 8^25) • (17.4) 

This demonstrates that both Lig and L25 are needed to arrive at a renormahza- 
tion group invariant formulation of the effective theory at next-to-next-to leading 
order. 



18 Matching U(3) and SU(3) 

The SU(3) framework applies if the quark masses and the external momenta are 
taken small compared to the mass generated by the topological susceptibility. 
Quantitatively, the condition takes the form W^ 



m, |(0|mm|0)| Otgd . (18.1) 

In this region, we may use the straightforward expansion in powers of momenta 
and quark masses also for the U(3) theory, so that the two effective descriptions 
have a common region of validity, on which we can compare them. 

As explicitly demonstrated in ref. the path integral for the effective U(3) 
Lagrangian indeed reduces to the one of SU(3) if the singlet meson field is inte- 
grated out. The calculation is analogous to the one described in [^], where the 
SU(3) effective theory is matched with the one relevant for SU(2). Like in that 
case, loops involving the propagation of light as well as heavy mesons require spe- 
cial attention, because the momentum scale of these is set by the heavy masses, 
while in the light sector, that scale does not occur. We briefly sketch the essential 



steps, referring to for a more detailed discussion. 

To match the two effective theories, we again use the decomposition intro- 
duced in section |1^: U = e^3^ U, ip = ip + 6. The phase factor explicitly exhibits 
the dependence of U on the singlet field ip and converts the relation det U = e^^ 
into the constraint det ?7 = e~*^. The mapping ensures that, under chiral ro- 
tations, U transforms in the same manner as U . Expressed in terms of these 
variables, the first term in takes the form 



33 



In the common domain of validity of the expansions in powers of p and 6, which 
is characterized by the inequahty ( p.8.1|) , the equation of motion imphes that the 



field ip represents a term of order p^, so that we may expand the expressions in 
powers of Tp. 

At leading order, the matching reduces to a comparison of the classical actions. 
Collecting the relevant pieces of and C^^\ we find that the effective U(3) 
Lagrangian reduces to 

C^Sf = \F'{D,WD^U) + + x^U) 

In the notation used here, the leading order SU(3) Lagrangian of ref. read^ 
C'^' = \F^D^U^D^U) + \F^{U^x + X^U) + ^(^o''' + F^)D^eD''e . (18.2) 

Hence the two theories match at leading order, provided the coupling constants 
F and B are the same in both versions and the couplings Ho are related by 

Hi""' = Ho + F^Ai + Oil/N,) . (18.3) 

In order to match the coupling constants of next-to-leading order, we need to 
perform the integration over the field ip, which describes the t]'. The key obser- 
vation here is that C^'^'^ contains the derivatives of this field exclusively through 
the term Y^F^{dip — DOy. This implies that - if the vacuum angle and the sin- 
glet axial field are turned off - only the vertices proportional to x generate loops 
involving the propagation of an rj'. Accordingly, the matching of the derivative 
terms is trivial: The coupling constants F,Li, ... , L^,L^,Liq and Hi are the 
same in the two versions of the theory. For B, Lg, -^7, Lg and if2, the matching 



In ref. [Q, the c ovariant derivative is defined as V^t/ — d^U + i {v^ + a^j)U — iU{Vf_i — a^). 



In view of eq. (12.3), this amounts to V^C/ = DfjU — ^Df^6U 
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conditions^ read 



5-3 = B|l-^Ao + 0(iV,-^)} 
^r-^a + ^(Ao-^)+0(iV-) 

3 = + L, + O(iV-i) (18.4) 



-SU3 



288 r 



Ao = A| , M2- 



6 r 



M-Mo ' 'J F2(l + Ai) 

The quantity Mq is the mass of the rj' in the chiral hmit. This mass sets the 
scale of the logarithm contained in Aq. As a check, we note that in the expression 
for Mo, the dependence of Ai on the running scale of QCD cancels against the 
one of r, in agreement with the fact that the mass of the rj' is renormalization 
group invariant, also in the chiral limit (this is the reason for not expanding the 
denominator in that expression). The difference between B^^'-^ and B accounts 
for the fact that the latter picks up a renormalization from the one loop graphs 
of the effective theory, while the former does not. 

The difference between the values of Lg, Lg, H2 in the two versions of the effec- 
tive theory is related to the difference between the corresponding coefficients Fg, 
Fg, A2. Indeed, one may check that in the matching relations between the renor- 
malized coupling constants, the divergences drop out. The coupling constant Lj 
is scale independent in both versions. In the extended theory, L7 is suppressed by 
the OZI-rule, while L|"3 represents a term of order N^. The leading contribution 
arises from T^'-exchange and is inversely proportional to the square of the mass of 
the i II : 

7-SU3 ^ 
-^7 - 



According to eq. (|18.4|) , the first order correction to this formula is determined 
by the OZI- violating coupling constants Ai and A2: 



^Note that the low energy constants are independent of the quark masses. The matching 
conditions do, therefore, not involve M^, Mk, Mjj, but they do contain the mass scale set by 
the topological susceptibility, which is related to the value of M^/ in the chiral limit. 
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The term Lj only enters as a correction of second order in Note that the 

dependence on the running scale of QCD also cancels out here: The relation 
( 16.7| ) shows that the ratio (1 + / (1 + Ai) is renormalization group invariant. 

In the present section we considered those couplings of next-to-leading order 
that survive if 9 and aP^ are turned off - the effective Lagrangian of ref. ^ only 
accounts for these. As shown in appendix |B|, the general expression of next- 
to- leading order involves 11 additional terms. The matching relations for the 
corresponding coupling constants are also given there. 



19 Anomalies 

Above, we focussed on the natural parity part of the effective Lagrangian. The 
photonic decays vr'^ 77, 77 77, rj' 77, for instance, are not covered, 
because these are contained in the unnatural parity part, which collects those 
terms that involve the tensor e^^p" , We now extend the above discussion to these 
and first consider the anomalies, which we analyze by means of the differential 
forms 

V = dx^v^ , a = dx^a^ , r = v + a , I = v — a , d = dx^d^ . 

The quantities dx^, dx^, dx"^, dx^ are treated as Grassmann variables. Their prod- 
uct yields the standard volume element, dx^dx'^dx'^'dx" = e'^'^^'^d'^. 

The phase of the determinant of the Dirac operator is not invariant under an 
infinitesimal U(3)j{XU(3)£ transformation of the external fields, 

Vr = 1 + z , = 1 + ^ ■ 

The part of the phase change that depends on the gluon field is compensated by 
the transformation 56 = (a^ — a^) of the vacuum angle. The remainder is unique 
only up to contact terms formed with the external fields v^, and 9. In the 
standard convention, it is invariant under the transformations generated by the 
vector charges, so that the change in the effective action is of the form 

55'e//{w,a, s,p,6'} = j - a^) Vt) . 
An explicit formula^ for Vt was given in ref. |^ : 



= ^ [F^F^ + \D^aD,a + \ i {F^ + AaF^a + a^F^) + \a^] , 
Fy = dv — iv"^ , D^a = da — iva — iav . (19-1) 



^The sign of fl is convention dependent; we use the metric H , set £0123 = +1 and 

identify 75 with 75 = -^70 7172 73- 
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The expression for Q is not unique, however - in fact, the one given is in conflict 
with the invariance of the effective action under the renormahzation group. The 
problem arises, because the singlet component of the axial external field trans- 
forms in a nontrivial manner under this group. The differential form a = dx'^a^ 



consists of two parts (compare section O) 



a = a + \D9. (19.2) 

The first is renormahzation group invariant and transforms as a gauge field un- 
der chiral rotations. The second transforms with [DOy^"^ = Z^^DO under the 
renormahzation group, but is invariant under U(3)iiXU(3)£. Inserting the de- 
composition in the expression for fi, we obtain a sum of three terms, 

VL = VLq + VLi + 0.2 ■ 
The first is obtained from Vl by replacing a with a: 

fio = ^ {KF^ + \D^aD^a + 1 1 {F^d? + AaF^a + a^F^) + \a^] , (19.3) 

where D^a = da — iva — iav. It is renormahzation group invariant. The remaining 
two terms are given by 



^i = J^{D.a{da)-i{F,a-dF,)De}, ^2 = {da) {da) , (19.4) 



and transform with fi^*^" = Z^^Qi, ^2^"' = Z^^Q2- If we were to identify the 
anomaly of the effective action with Q, this functional would fail to be independent 
of the running scale of QCD, in contradiction with eq. (|15.1|) . 



The problem is readily solved: The extra terms Qi and Q2 represent the 
anomalies generated by two contact terms, 

Pi = ,^{aD,a)D9 , P, = {a){da)d9 . (19.5) 

Removing these, the anomaly takes the renormahzation group invariant form 

5Seff{v,a,s,p,9} = j{{aL - ««) ^0) • (19.6) 

Since the contact terms Pi, P2 vanish for (a^) = df^9 = 0, they only matter when 
considering correlation functions that contain the operators ^4°, u. 
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20 Wess-Zumino-Witten term 



Within the effective theory, the anomahes are accounted for by the Wess-Zumino- 
Witten term. The standard expression for this term reads 



W{U,rJ) = {Ul^U^r + \UlU^rUlU^r + iUdllU^r + idrUlU^ 



-i^lU^rUl + ^ UUr U I - S^t/V Ul + ^ldl + ^dll 

+ \i:iT.l-iY:H) -{R^ L) , (20.1) 

with U G U(3) and S = WdU. The first term is an integral over a field U{x, x^) 
that smoothly interpolates between U{x,0) = 1 and U{x, 1) = U{x). The Grass- 
mann algebra is supplemented with a fifth element dx^ and the integration ex- 
tends over the five dimensional manifold M5, which represents the direct product 
of Minkowski space with the interval < < 1. The integral is independent 
of the particular interpolation chosen to connect U{x) with the unit matrix. In 
the second term, the integration only extends over Minkowski space, M4. The 
operation {R ^ L) requires an interchange of the 1-forms r and / as well as an 
interchange of U and U'^ . By construction, an infinitesimal chiral rotation of the 
variables [/, v and a generates the change 

5^wzw{f/, v,a}= [ {{a, - a^) n) , (20.2) 



where Q is the 4- form specified in eq. ( |19.1|) . 

For the reason given in the preceding section, we remove the contact terms Pi, 
P2 and define the Wess-Zumino-Witten part of the effective Lagrangian through 

jdxC^^^ = 5wzw{t/, v,a}- j (Pi + P2) . (20.3) 

This ensures that the term £wzw accounts for the anomalies of QCD in the 
renormalization group invariant form ( p.9.6D : 



6jdxCy,^y, = j {{a r,- a ^) no) . (20.4) 

Note, however, that the term £wzw as such is not renormalization group invari- 
ant. Since the scale dependent pieces contained therein are gauge invariant (see 
appendix 0), we could remove these and arrive at a scale invariant expression. 
The drawback of such a choice is that it interferes with the large counting 
rules: As demonstrated in the next section, our definition of C^zw is singled 
out by the property that it represents the leading unnatural parity part of the 
effective Lagrangian. 
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21 Unnatural parity part beyond leading order 



We denote the unnatural parity part of the effective Lagrangian by £e//- As we 
did not find a complete list for the terms of order in the literature, we briefly 
outline the construction. 

Once the WZW-term is removed, the unnatural parity part also becomes 
gauge invariant under 11(3)^x11(3)^,. It is convenient to express the Lagrangian 
in terms of the variables U, ip, v^, a^, D^d and their derivatives. Terms involving 
derivatives of if) can be integrated by parts. At order p^, charge conjugation 
invariance then allows six independent invariants i^F^^ = ^e'^"'"^ Fp^j): 

+v, {{Rn{R,u) + {Ln{L,u)) + % {Rn{L,u) + o{p') . 

On account of parity, all of the potentials are odd functions of tp, except for V4, 
which is even. 

In the 1/Nc expansion, the leading contribution to the potentials Vi, V2 and 
V3 is linear in ip and contains a coupling constant of 0(1), while V4 reduces 
to a constant that is also of this order. Since the remaining two terms involve 
two flavour traces, their expansion starts at order This implies that the 

simultaneous expansion in powers of p and is dominated by the WZW-term, 
which represents a contributions of order NcP^ = 0{6). The remainder is of order 
5"^ or higher: 

4// = >Cwzw + + £(3) + . . . (21.1) 

The general expression for the next-to-leading order Lagrangian contains con- 
tributions of the type and NcP^. The former can be extracted from the repre- 
sentation given above: 

C^^ = L, %^ {R^'^DpUD.U^ + Lf^'-D^U^DM) + {R^'UL^.U^) (21.2) 
+L3 {R'^'R^.. + ^'^^ V) + ^4 iD^e {Rf^'^D^UU^ - L^'^U^DM) . 

At order N^p^, many invariants can be formed, in particular also terms propor- 
tional to the quark mass matrix. Below we will discuss only a selection thereof: 
The terms relevant in connection with the radiative transitions. 

Note that £wzw is not independent of the scale used to renormalize the sin- 
glet axial current. The scale dependent part is gauge invariant, but represents 
a contribution of leading order in the 1/A^ expansion and must be retained for 
the relation ( |21.1|) to hold. As in the case of the natural parity part, renormal- 



ization group invariance thus requires specific contributions of nonleading order. 
The phenomenon arises from the fact that the loop graphs responsible for the 
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anomalous dimension of the singlet current violate the OZI rule, which implies 
that some of the effective coupling constants contained in must compensate 
for the scale dependence of i2wzw The relevant terms are those of order p^, all 
of which were listed above. The dependence of the coupling constants on the 
running scale of QCD follows from the decomposition of the WZW-term given in 
appendix 0: The sum £wzw + '^p4^ is renormalization group invariant to order 5^, 
provided Li, . . . , L4 are renormalized according to 

L4 =Z^L^ + K, K = —-^-^. (21.3) 
The quantity n is of 0(1), like the coupling constants themselves. 



22 Radiative transitions 

As an illustration, we consider the radiative transitions vr^ — > 77, r] — > 77, 77' 
77 |3^. At leading order, the corresponding part of the effective Lagrangian is 
obtained from the Wess-Zumino-Witten term with 

r = 1 = -eQA , 

where Q = diag{|, — |, — i} represents the charge matrix of the light quarks 
and A = dx^A^ is the 1-form associated with the electromagnetic field. Setting 
[/ = e**^, the terms linear in and quadratic in A reduce to 

^wzw{f/,t^,a} = ^ j{Q^d<p)AdA=-^ jd'x{Q'<P)F,,F^'' . 

Concerning -^^4 , only the combination L2 + 2L3 of coupling constants matters 
for the photonic transition matrix elements. The net effect of the first order OZI 
violations is that the trace {Q'^cj)) appearing in the WZW-term is replaced by 

withQ Ki = — 487r^(L2 + 2L^)/Nc. The scaling laws for the coupling constants L2, 
L3 imply that the renormalization of Ki is of the same form as the one for A2: 

i+i^r" = ^^(i+iri) . (22.1) 

One readily checks that this indeed compensates the renormalization of the singlet 
field in eq. ( |16.2|) , so that the result for the transition amplitude is renormalization 
group invariant. 



^in ref. Ilq] this coupling constant is denoted by A3 = Ki. 
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The Lagrangian C'^'^'^ contains two further categories of contributions: Terms 
of O(A^) with 6 derivatives and chiral symmetry breaking effects of O(A^mp^). 
We denote these by C-^J and , respectively. As far as the photonic transitions 
are concerned, the first contains two independent contributions, which may be 
written in the form: 

Using the equations of motion, the second term can be absorbed in the couphngs 
occurring in £^4 and Cx ■ Since the term with DF^*^ only matters for off-shell 
photons, we can ignore this part of the Lagrangian altogether. 

Finally, we consider the chiral symmetry breaking terms. As only neutral 
mesons are involved, the matrices 0, Q and x = 25m commute, so that there is 
only one independent invariant with a single flavour trace. We again extract a 
normalization factor and denote the coupling constant by K2: 

The net result for the effective Lagrangian that describes the photonic decays to 
first nonleading order thus reads: 

Cp^,, = {{Q'<P) + \Kim{<t>) + MQ\<P)] F,,F^^ . (22.2) 

The coupling constant Ki describes the corrections generated by the violations 
of the OZI rule and K2 accounts for the breaking of chiral symmetry. While 
Ki represents a term of order and depends on the running scale of QCD 
according to eq. (|22.1|) , the constant K2 is of order 1 and is renormalization group 
invariant. 



23 Summary and conclusion 

The effective theory of QCD with three colours is well known. In that framework, 
the effective Lagrangian consists of all terms respecting chiral symmetry. The 
chiral perturbation series amounts to an expansion in powers of momenta and 
quark masses. In the present paper, we have examined the extension of this 
effective theory required to analyze the low energy properties of QCD in the limit 
where the number of colours, N^, is treated as large. In that case, the situation is 
more intricate, because there is an additional low energy scale, related to the mass 
of the 1]'. As discussed in detail, the standard expansion in powers of momenta 
and quark masses, which we refer to as the p-expansion, cannot be used when 
becomes large: In order to coherently analyze the low energy properties of the 
T]', the two quantities appearing in the r^'-propagator - the mass and the square 
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of the momentum - must be treated on equal footing. We exploit the fact that 
M^, represents a term of 0(1/A^), so that an expansion that counts p^, m and 
1/Nc as quantities of the same order does have the desired property. We refer to 
the corresponding low energy expansion as the 5-expansion. It orders the triple 
series in p = O(v^), m = 0{5) and = 0{5) by collecting terms that are of 
the same order in 5. 

In the limit Nc —>■ oo, g'^Nc fixed, the perturbative analysis of QCD determines 
the order in 1/Nc of the various correlation functions, according to eq. ( p.3| ). 
These counting rules only hold for nonexceptional momenta. In particular, the 
low energy singularities generated by the Goldstone bosons (including the r]') give 
rise to contributions that violate these rules. We have shown, however, that the 
generating functional Seff{v, a, s,p, 6} does admit a coherent expansion in powers 
of (5. 

What we then discussed in some detail are the implications of this property of 
the generating functional for the effective Lagrangian. The reason why this step 
is not straightforward is closely related to the fact that the dynamical variables 
of the effective theory do not have physical meaning - their choice is inherently 
ambiguous. This entails that statements about the effective Lagrangian can only 
be true modulo a change in the dynamical variables. What we were able to show 
is that there exists a class of such variables, for which the coupling constants 
in the effective Lagrangian are at most of 0{Nc). This is natural, because the 
correlation functions collected in the generating functional obey the same bound. 
Furthermore, we have shown that, like for the generating functional, the depen- 
dence on the vacuum angle 6 is suppressed in the large limit: At order 5", the 
effective Lagrangian is a polynomial in 6, of degree n + 2. 

The outcome of our investigation boils down to a remarkably simple construc- 
tion recipe for the effective Lagrangian that holds to any given order in 6: 

(i) Apart from the WZW-term, the Lagrangian is manifestly invariant under 
local U(3)fl.xU(3)i transformations. 

(ii) At order 6", the Lagrangian represents a polynomial formed with the fields 
U, W , s, p, R^^, L^^, ip, D^6 and their covariant derivatives. The expression 
contains terms that are at most of order p2n+2_ 

(ill) The coupling constant associated with a term in the Lagrangian is of order 
N^~^ , where k counts the number of traces plus the number of factors with 
D^9 or derivatives thereof (note that k>l). 

In the text, we explicitly give the expressions for the Lagrangians of order 1 
and order 5, but list only a selection of the vertices occurring at next-to-next-to 
leading order, for practical reasons: We expect the full Lagrangian of order 5^ to 
contain about 100 terms. 
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We have also shown that the effective theory is consistent with periodicity in 
the vacuum angle 6. This is not evident, a priori, because the effective Lagrangian 
truncated at a given order in 6 is actually a polynomial in 9. The paradox 
is resolved by observing that the periodicity is not a necessary feature of the 
Lagrangian itself - only the corresponding effective action needs to be periodic, 
and this is the case. 

A corollary of our results concerning the structure of the effective Lagrangian 
is that the transformation of Kaplan and Manohar is forbidden at large A^: For 
6^0, this transformation involves the vacuum angle through a factor e~*^ and 
thus generates a modification of the effective Lagrangian with a nonpolynomial 
dependence on 9. This is in confiict with the properties of the effective theory at 
large N^. 

As is well known, the singlet axial current A° carries anomalous dimension 
and thus depends on the running scale of QCD. Moreover, the renormalization 
group mixes the operators uj and In the effective action, the correspond- 

ing external fields are the trace (a^) and the vacuum angle 9. We have shown 
that their renormalization group properties follow from symmetry considerations 
alone: While the vacuum angle is scale independent, the singlet field (a^) trans- 
forms inhomogeneously under scale transformations and picks up a contribution 
proportional to the gradient of the vacuum angle. We have worked out the con- 
sequences for the dynamical variables and coupling constants of the effective 
theory. In the leading order Lagrangian the coupling constant F is scale 
independent, while B and r are multiplicatively renormalized. Some of the fields 
contained therein, however, transform inhomogeneously, so that C^^^ does not re- 
main invariant when the running scale is varied. The origin of the problem - the 
anomalous dimension of the singlet axial current - is due to graphs of nonleading 
order. This implies that the change in produced by a change of scale is an 
effect of order 1/Nc and is eaten up by a suitable shift of the couplings occurring 
in C^^\ so that the effective action does remain invariant. Quite generally, the 
action of the renormalization group on the effective coupling constants occurring 
at a given order mixes these with the lower order couplings. Physical quanti- 
ties only involve scale independent combinations of coupling constants, which are 
easily identified. 

In the remainder of the paper, we have extended the analysis to the unnatural 
parity part of the Lagrangian and have shown that the scale independence of the 
effective theory can be made manifest also in this sector. At leading order, the 
unnatural parity part is given by the Wess-Zumino-Witten term, which accounts 
for the anomalies within the effective theory. The straightforward extension of 
this term to the case of U(3)ijXU(3)£,, however, fails to be invariant under the 
renormalization group. In part, the deficiency only concerns contact contributions 
that are readily removed. The remainder still contains scale dependent contri- 
butions, but these are gauge invariant, so that the modification produced by a 
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change of the running scale may be absorbed in the couphng constants occurring 
at the next order of the expansion. The relevant terms are those of unnatural 
parity at order 5^ . We could instead have modified the expression for Cwzw, by 
adding suitable gauge invariant pieces that make it scale independent. This, how- 
ever, would upset the large Nc counting rules, so that the leading contribution in 
the 5-expansion of the unnatural parity Lagrangian would then not coincide with 

■^wzw- 

To establish contact with the standard low energy theory of QCD, we have 
given the explicit matching relations between the low energy constants relevant at 
large and at fixed A^. We have also discussed the extension needed to investigate 
singlet currents in the standard framework. Some features, such as the depen- 
dence of the effective coupling constants on the running scale of QCD apply to 
both versions of the theory. Others do not: The Kaplan-Manohar transforma- 
tion, for instance, is in conflict with the properties of QCD only if the number 
of colours is treated as large - at fixed the effective theory is invariant under 
this operation. 

A Construction of the effective Lagrangian 

In the present appendix, we show how the reasoning of sections can be ex- 
tended to construct the full effective Lagrangian. As input, we use the large 
properties of the Green functions of QCD. We first switch the quark masses off 
and use the large counting rules for the scattering amplitudes of the pseu- 
doscalar mesons. These allow us to establish corresponding counting rules for the 
interaction vertices of the effective Lagrangian. Then, we generalize the argument 
to those vertices that describe the response of the system to the perturbations 
generated by the external fields and finally discuss the consequences of the Ward 
identities of chiral symmetry. 

Scattering amplitude at large 

If the number of colours is sent to infinity and the quark masses are turned off, 
the spectrum of QCD contains nine massless pseudoscalar mesons. The scattering 
amplitudes describing the interaction of any number of these particles in the 
initial and final state can be extracted from the connected correlation functions 
formed with the corresponding number of axial currents, using these currents as 
interpolating fields. The correlation functions represent quantities of order A^, 
irrespective of the number of currents contained therein. In view of the fact that 
the one particle matrix elements of the currents, (0|y4^|'7r''), are of order y/N^, the 

scattering amplitude for n = rii + Uf mesons is at most of order Nc~"'^'^. 

The scattering amplitude contains singularities in the low energy region. In 
particular, for n > 6, it contains one particle reducible contributions, describing 
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a sequence of collisions, connected by the exchange of single particles. Denoting 
the number of exchanged particles by there are i + 1 irreducible parts. The 
number of meson lines entering or leaving the irreducible parts adds up to n + 2i. 
Applying the counting rule to the irreducible parts, the resulting contribution 
to the scattering amplitude represents a term of order N^, with k = (£ + 1) — 
I (n + 2i) = 1 — i ra. This shows that the singularities generated by one particle 
exchange manifest themselves already at leading order. 

Unitarity relates the imaginary part to the square of the scattering ampli- 
tude. The relation implies that this amplitude contains further singularities. The 
exchange of a pair of mesons (two particle intermediate states in the unitarity re- 
lation), for instance, generates a branch cut in the one-particle-irreducible parts. 
The contribution from the two particle cut, however, only shows up at order 
Nc : If the two parts connected by the two exchanged particles involve rii and 
n2 mesons, respectively, we have rii + n2 = n + 4, so that the overall power of 

is given by (1 — ^) + (1 — ^) = — |. Exchanges of more than two parti- 
cles between the same irreducible parts are suppressed even more strongly. This 
means that the branch cuts required by unitarity only show up at nonleading 
orders of the expansion. The leading order contributions only contain those 
singularities that arise from one-particle-exchange. Moreover, at leading order, 
the one-particle-irreducible parts reduce to polynomials of the momenta. 

At leading order of the 1/Nc expansion, the structure of the scattering ampli- 
tude is the same as the one of the tree graphs of a pseudoscalar field theory. We 
identify the dynamical variables with the dimensionless fields . . . , 

introduced in eq. ( p.3|) and represent the interaction Lagrangian in the symbolic 
form 

Ceff = J29ik,n)xd'xr , (A.l) 

k, n 

where the flavour and Lorentz structure of the vertices is suppressed. The integers 
k and n merely count the number of derivatives and fields occurring in the vertex 
in question and g{k, n) represents the corresponding effective coupling constant - 
in general, there are several, independent vertices of the same symbolic structure. 
Lorentz invariance implies that k is even. 

The terms quadratic in 0(x) are given in eqs. ( |6.2|) and ( ^.4|) . The correspond- 
ing coupling constants are 5'(0,2) ~ = 0(1) and 5'(2,2) ~ = O(A^). For 
n > 2, the coupling constant g{k,n) generates a tree graph contribution to the 
one particle irreducible scattering amplitude with n mesons. The contribution 
is of the symbolic form g{k,n) p'^ F~^, where p stands for the momenta of the 
particles. The comparison with the counting rule for the scattering amplitude 
suggests that the coupling constant g{k, n) can at most be of order A^. 
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Freedom in the choice of the dynamical variables 



Actually, the argument just given runs in the wrong direction: It only shows that 
if the Lagrangian exclusively contains vertices of order A^, then the corresponding 
scattering amplitude does obey the large counting rule - the converse is not 
true. A counter example can be constructed as follows. As discussed in section 
1^, the effective Lagrangian is not unique, because its form depends on the choice 
of variables. We may for instance subject the field U{x) to the transformation 
U' = U expiflip + 6). The operation preserves the transformation law 1| ) , 
irrespective of the choice of the function f{x). We may choose one that grows 
with A^. Suppose that the coupling constants are of order A^ and express the 
Lagrangian in terms of the new variables. The resulting expression describes the 
same physics, but contains effective coupling constants that grow more rapidly 
than with the first power if A^ becomes large. 

The freedom in the choice of the dynamical variables is related to the fact 
that, in the scattering amplitude, all of the momenta are on the mass shell. It 
is well known that different vertices may give rise to the same on-shell matrix 
elements. The requirement that the on-shell matrix elements of a given set of 
vertices reproduces certain contributions occurring in the scattering amplitude 
only fixes these matrix elements up to terms that vanish on the mass shell of 
the colliding particles, = Mf. In the chiral limit, eight of these are massless. 
Mi = 0, while the mass of the ninth is given by M^/. The off-shell extension 
involves an ambiguity of the form YliAvf ~ ^f) Ci{Pi, ■ ■ ■ Vn)- In coordinate space, 
this ambiguity corresponds to terms in the Lagrangian that are proportional to 
the equation of motion, which at leading order in the 5-expansion of the massless 
theory is of the form □ + f((/)) = with f = 2r/F^. The right hand side, 

consists of a series of terms that contain three or more fields. 

This observation may be used to determine the ambiguity in the effective 
Lagrangian in an iterative manner. Suppose that the tree graphs of Ceff and C^j^ 
generate the same on-shell scattering matrix elements. As discussed in section 
1^, we may choose the variables such that the terms that are quadratic in the 
meson fields are the same, so that the difference A£ = C'^^^ — Ceff only contains 
vertices with four or more meson fields. Suppose now that the four-particle 
scattering amplitudes generated by these Lagrangians coincide. This property 
implies that, up to a total derivative, the terms of order cf/^ contained in A£ 
can be written in the form c{k) x (□(/) + r(0)) x 5^ x (p^. We may replace 
□ + f(0) with □(/) + f (0) — because the extra contributions contain six 
or more meson fields. It therefore suffices to transform the variables in C^^j with 
(j) + F~'^ c{k) X 9^ X 0^: At order 0^, the operation reduces the difference 
between the two Lagrangians to a total derivative. Since a change of variables does 
not modify the physics, the new version of C^^j yields the same scattering matrix 
elements as the original one, irrespective of the number of particles participating 
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in the collision. Iterating the procedure, we may extend this analysis to terms 
with an arbitrary number of fields. We conclude that the on-shell scattering 
matrix elements unambiguously determine the effective Lagrangian, except for 
two inherent degrees of freedom: Choice of the dynamical variables and total 
derivatives. The first refiects the ambiguities occurring in the extension off the 
mass shell, the second concerns the extension off the energy-momentum shell. 

We will fix the choice of variables when specifying the explicit expressions for 
the first few terms of the derivative expansion. For the moment, we only exploit 
the fact that there is a set of coupling constants g{k,n) = 0{Nc), for which the 
tree graphs of the effective Lagrangian do reproduce the scattering amphtudes at 
leading order of the l/N^ expansion. This only excludes those transformations of 
variables that generate coupling constants growing more rapidly than with the 
first power of A^. 

External fields 

Let us now turn on the external fields. The above analysis is readily extended to 
this case. We again consider the massless theory and use the symbol ji = q Tiq 
to denote any one of the quark currents. The large counting rules for the 
correlation functions of these operators were given in section ^. We first generalize 
this to matrix elements between asymptotic states and consider 

Gn.n^n^Uf = {f \ T ji{xi) ■ ■ ■ jn^^Xn^) Uj{jji) ■ ■ ■ Uj{yn^)\ l) c , 

where \i) and (/| represent states with incoming and nj outgoing mesons, 
respectively. The behaviour of Gn^n^n^uf at large is established in the same 
manner as for the scattering matrix: The matrix element is related to the residue 
of the poles occurring in a correlation function that, in addition to the operators 
listed, involves n = rii + nj axial currents, which play the role of the interpolating 
fields, while the operators ii{xi) ■ ■ ■a;(?/„^) are treated as spectators. Denoting 
the contribution to Gn^n^mnj that arises from graphs with I quark loops by 
Gil- m rif 5 the generalization of eq. ( ^TSf ) reads 

Gl,..n..,=0(iV^'-"--^'^^-^"^) , £=1,2,... (A.2) 

We may also check that, like in the case of the scattering amplitude, the singu- 
larities generated by one particle exchange manifest themselves at leading order, 
while the unitarity cuts only appear at nonleading orders. In the large Nc limit, 
the one particle irreducible parts can therefore again be expanded in the mo- 
menta. 

The above counting rule shows that the matrix elements of all of the quark 
currents behave in the same manner in the large limit, while those of the 
operator u{x) are suppressed. For the counting of powers relevant at low en- 
ergies, on the other hand, the external field 6{x) counts as a term of order 1, 
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'x),a^(x) = 0{-\f5) and s(x),p(x) = 0{5). In the following bookkeeping, we 
do not distinguish between v^{x) and a^(x), nor between s{x) and p{x) and write 
the effective Lagrangian in the symbolic form 

C^^j = ^e(A;, ke) x N^^~'"> x d'' x (A.3) 

+ ^' '^^'f ^ ^c'''' X '9'' X 0" X t;'^" X s''^ X e'^^ . 

The first sum accounts for the contributions generated by graphs that do not 
contain quark lines. It represents the effective Lagrangian of gluodynamics, which 
we discussed in section |[ The integers k and kg count the number of derivatives 
and the number of times the external field 9{x) occurs, respectively. In eq. (p.7[), 
that Lagrangian is written in the form —N^eo{'&) + d^dd ^dei{d) + . . . , with 
d = 6/Nc. As we are now expanding in powers of the vacuum energy 

density N^eo{{}) is replaced by the series A^^eo(O) + | 6*^60(0) + . . . and likewise 
for the other coefficients. The comparison shows that the coupling constant 
e{k, kg) is at most of order 1. On account of parity, e{k, kg) vanishes unless 
kg is even. Disregarding the term e(0, 0) = — A^^eo(O), which merely contributes 
to the cosmo logical constant, the sum only starts at A;^ = 2. 

The second part of the effective Lagrangian arises from graphs containing at 
least one quark loop. The integers k and n count the number of derivatives and 
meson fields, k^ is the number of external vector and axial fields, kg counts the 
scalar and pseudoscalar ones, while kg is the number of times the field 6{x) enters. 
As discussed above in connection with the scattering matrix, the translation of 
the counting rule for the matrix elements into one for the vertices of the effective 
Lagrangian involves ambiguities related to the freedom in the choice of the dy- 
namical variables. In the presence of external fields, that freedom becomes even 
richer, because the variables may be subject to transformations that depend 
on these fields. It suffices to observe, however, that the one particle irreducible 
matrix elements are polynomials in the momenta, so that their Fourier transforms 
represent a collection of delta-functions and derivatives thereof. We may simply 
multiply this object with the relevant external fields and add a factor of 
for each one of the on-shell mesons. Integrating all but one of the coordinates over 
space, we obtain a specific representation for a term in the effective Lagrangian, 
for which the relevant tree graph does reproduce the matrix element in question. 
The effective coupling constants g{k, n, ky, kg, kg) occurring therein are the co- 
efficients of the polynomial that describes the matrix element. Their order in the 
1/Nc expansion is determined by the counting rule (|A.2|) , which thus implies that. 



if the effective Lagrangian is constructed in this manner, the coupling constants 
g{k, n, k^, kg, kg) are at most of order 1. The Lagrangian in eq. (|A.1| ) is what 
remains if all external fields are turned off: g{k, n) = Ncg{k, n, 0, 0, 0). 
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Chiral symmetry 

We now consider the simultaneous expansion in powers of momenta (or deriva- 
tives) and introduced in section The above symbohc expression for the 
effective Lagrangian exphcitly displays the number of derivatives, but only indi- 
cates the leading power of A^: The expansion of the effective coupling constants 
in powers of starts at 0(1), but also contains terms of nonleading order. 
In the following, we stick to this abbreviated notation. With the assignments 
specified in eqs. ( |10.1| ), ( |10.2| ), the general vertex represents a term of order 

Lorentz invariance implies that k + is even, so that only integer powers of 6 
occur. Ordering the Lagrangian in this manner, it takes the form 

where the term collects all contributions of 0(5"). Actually, as shown below, 
chiral symmetry implies that the first term vanishes. The expansion only begins 
with the term C^'^\ which collects the contributions of 0(1). Note that only a 
finite number of derivatives and external fields can occur at any finite order in 6, so 
that the relevant Lagrangian only involves a finite number of coupling constants, 
like in the standard framework. In particular, at 0(5"), the dependence on the 
vacuum angle is a polynomial of degree 71 + 2. This feature reflects the fact that, 
in the large limit, the 9 dependence is suppressed, both in gluodynamics and 
in massless QCD. 

As mentioned in section 0, the Ward identities of chiral symmetry are equiva- 
lent to the statement that the effective action is invariant under the 11(3)^1 xU(3)i 



gauge transformation of the external fields specified in eq. (2^). Note that the 
transformation law only relates quantities of the same order in 6. The effective 
Lagrangian can therefore be gauge invariant only if this is the case separately 
for each one of the terms In particular, the contributions of order 

must altogether represent a gauge invariant expression. Now, these arise from 
k = ky = kg = kg = and thus only depend on the meson field. We may think 
of this part of the effective Lagrangian as being a function of the field U{x), 
which, moreover, does not involve derivatives thereof, = f{U). Invariance 

under U(3)«xU(3)i implies that this function obeys fiynUVl) = f{U). With 
Vj, = U .,Vr = 1, this leads to f{U) = const. Hence C^~^^ only contributes to the 
cosmological constant and may be discarded. 

It is convenient to replace the vacuum angle hj ip = ip + O and to use the fields 
U,v,a,s,p,'ilj and their derivatives as independent variables - the Lagrangian £(") 
is a gauge invariant function thereof. Moreover, this function is a polynomial in 
all variables except U. Note that only the fields at one and the same point of 
space-time enter, so that it is legitimate to treat the derivatives df^^dyU . . . 
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as independent from U. Since ip by itself is invariant under U(3)jixU(3)i, gauge 
invariance does not restrict the dependence on this variable and its derivatives. 

At leading order, the above counting rule permits five independent invariants^: 
{WDf'D^U), {D^UWf'U), {{s + ip) W), {{s - ip) U), V-^. The first term differs 
from the second only by a total derivative and can thus be discarded. For the 
Lagrangian to be real, the coefficients of {{s + ip)W) and {{s — ip)U) must 
be complex conjugates of one another and parity then implies that they are 
real. Hence the leading order Lagrangian contains three independent coupling 
constants. The explicit expression is given in eq. ( |10.5| ). 

It is straightforward to generalize this procedure to find the expressions for 
the higher order Lagrangians. At order 5, for instance, the relevant terms are 
those of order p^ with two flavour traces and those of order with one trace. 
Using the equations of motion associated with C^^\ the Lagrangian of order 5 
may be brought to the form given in eq. ( |10.6D . 



B Full second order SU(3) Lagrangian 

In ref. 0], the effective Lagrangian relevant for the low energy analysis of matrix 
elements involving the winding number density uj{x) or the singlet vector and axial 
currents V^(x), was given only at leading order. In the present appendix we 

briefly discuss the extension needed to study these quantities to first nonleading 
order of the p-expansion. In the language of the effective theory, the relevant 
extension is obtained by introducing three additional external fields, 6* (a;), {v^{x)), 
(a^(x)) and replacing the symmetry group SU(3)hxSU(3)i, by 11(3)^x11(3)^. 

The leading order SU(3) Lagrangian is given in eq. ( |18.2| ). At order p"^, the 
terms listed in ref. [Q can be taken over as they are, simply replacing the variables 



U, and V^U with the quantities U, ttfj, and -D^?7, defined in section |T2 



/-sua 



+Lf%U\ - xW + LT'{U\U\ + x^Ux^U) (B.l) 

where R^u, L^y are the field strengths of ± a^. 

In addition, chiral symmetry, parity and charge conjugation invariance permit 
11 new couplings, containing the covariant derivative of the vacuum angle or the 



"^Note that the term {W D^^U) {W Dfj_U) involves two flavour traces and thus only occurs at 
first nonleading order of the expansion. 
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field strengths of tlie singlet external field^: 

-tLlf D^eiDi^Wx - D^Ux^) + tLlfd.D^e {Wx - X^U) 
+iL'^'t^'""'D^e{R,pD„UU^ - L.pU^D^U) (B.2) 
+Hp{R,, + + in + Hf^R,, - L,,){R^'' - in 

Finally, the anomalies of the underlying theory require an extra term that is 
not gauge invariant, but does not involve any free constants. In the notation 
introduced in eq. ( p0.1| ), this term is given by 



j dxC^^l^ = S^zw{U, V, a} . (B.3) 
The full effective Lagrangian of order reads 

c;? = cT + ^b' + ^wlw • (B.4) 

Formally, the extension of the low energy analysis to the matrix elements of 
the operators uj{x), A^^{x), V^(x) thus nearly doubles the number of effective 
coupling constants. Four of these, however, represent contact terms and a fifth 
only concerns matrix elements of unnatural parity (note that a coupling of this 
type does not occur in the Lagrangian of ref. - Lf^^ only matters for matrix 
elements that involve the winding number density or the singlet axial current). 

The inclusion of the scale dependent field (a^) implies that some of the cou- 
pling constants must be renormalized also in this version of the theory. The renor- 
malization procedure does, however, not entangle coupling constants occurring at 
different orders in the expansion. In fact, in the above basis, the renormalization 
is homogeneous: In contrast to a^, the field is renormalization group invariant. 
The constants Lf^^ , . . . , L\q^ are therefore scale independent, while the remain- 
ing couplings pick up a multiplicative renormalization that compensates the one 
of the external fields, 

I Tiren , t ren\ / p < t \ / Tjren j ren\ T \ 

The renormalizations needed to absorb the infinities generated by the one loop 
graphs of £^2^ may be worked out as follows. The change of variables IJ = es^U 

*In ref. | |3^ , only the coupling constants relevant for the two-point-functions were considered. 
For a corresponding list of U(3) invariants, see |14|. 
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takes that Lagrangian into 

det U = 1 , D^U = djj -i{v^, + d^)U + iU{v^- a^) , 

where and are the traceless parts of the external fields, Xg = ei^x ^"^^ 
Hf) = Hq^'^ + F^. In this form, the angle 6 and the singlet axial field exclusively 
occur in Xg and -D^^, so that the one loop calculation of ref. can be taken 
over as it is, simply replacing x by Xg. In particular, that calculation shows that 
divergences proportional to the new terms do not occur, so that the corresponding 
renormalization coefficients Fn, . . . , Tu and A3, ... , Ag all vanish. 

The matching relations for the standard part C^^^ of the SU(3) Lagrangian 
are given in section |T^. These specify the leading terms in the expansion of 
the coupling constants Lf^^, . . . , L^q^. For those in the part involving the singlet 
external fields, the analogous relations read[| 

Lff = -4(L2 + iL3) + 0(l) 
Llf = l{L, + lL, + lL,) + Oil) 
Llf = liL, + lLs)+Oil) 
Llf = liL, + 3L, + L,s) + 0il) 
Llf = l{2L, + 3L,s)+Oil) 

Llf = -F\l + A0(1 + A2) (72 r)-i + 0(1) (B.5) 
L\f = iVe(288 7r2)-i + \U + 0{N-^) 
HP = 0(1) 

i/r = |(ifi-iLio)+0(l) 
i/r = |(Li + iL2 + |L3) + 0(l) 
^su3 ^ ^4(^ ^ ^^)2 r)-i + 0(1) . 

The coupling constants Lfg^ and H^^^ receive a contribution from r^'-exchange, 
similar to the one in Lf^^ (section |18D. The leading contribution to Lf"^ stems 
from the Wess-Zumino-Witten term of the extended theory (appendix 0). 

Note that, in the SU(3) framework, the Kaplan-Manohar transformation ( [13. 1| ) 
takes the mass term of the effective Lagrangian into 

{u^x) - {u^x) + ^ {{xW - (x^ux^u)} , 

even if the vacuum angle does not vanish: The factor e~*^ is absorbed in the 
field U, on account of det U = e~*^. Hence the transformation of the quark mass 
matrix in eq. (|13.1|) is equivalent to a change of the effective coupling constants 



occurring in C^^, also for 6* 7^ 0. 

^Note that the independence with respect to changes in the chiral renormahzation scale is 
not manifest in these equations - the dependence of the couphng constants on the scale fi is an 



effect of order and thus beyond the given accuracy. 
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C Renormalization of the WZW term 



We first observe that the field U introduced in eq. (|12.2|) is renormahzation group 



invariant. Moreover, U and U transform in the same way under chiral rotations, 
because the factor that makes the difference, e^'^, is invariant. This immediately 
implies that the anomalies of the functional S'wzw remain the same if U is replaced 
by U. Indeed, performing the above change of variables in eq. ( p0.1| ), we obtain 
a gauge invariant result for the difference: 



Syvzw{U,v,a} = Swzw{U,v,a} + J A 



3 
UAtt 

F, = dr- ir^ , F, = dl - il^ . 



A = -77^ V {iFr DUDU^ + Fr UFi + 2F^ + {R^ L)) 



with DU = dU — irU + iUl. The terms occurring here are of the same form as 
those entering the U(3)flXU(3)£ invariant part of the effective Lagrangian. 

Next, we decompose the axial field according to eq. (|19.2|) . The corresponding 



decomposition of the right- and lefthanded gauge fields into a renormalization 
group invariant part f, T and a remainder reads r = f + \D9, 1 = 1 — \D9. Using 
the identity {dUU^) = —id6, which follows from det U = e~*^, we then obtain 

5wzw{^, V, a} = 5wzw{f/, v,a} + y (A + 5 + Pi + P2) , 
UAtt 

Fr = df-if^, Fi = dl-il^ 



B = iDO {Fr DUU^ - Fi U^DU) 



Note that terms proportional to DO {{DUV^Y) cancel out on account of charge 
conjugation invariance. The term S'wzwjt^; a} is renormalization group invari- 
ant. Under chiral rotations, it transforms with 



^^wzwlf^, v,a} = J {{a^ - aji)Qo) . 

The one with B is gauge invariant, but transforms with under the renor- 
malization group. The calculation automatically yields the contact terms Pi, P2 



introduced in section |19|, which account for the difference between f2 and f2o - 
these transform in a nontrivial manner, both under the renormalization group 
and under chiral rotations. Finally, A may also be sorted out according to the 
behaviour under the renormalization group: 

^ = -T77-^ ^ i^^r DUDW + Fr UFi + + ^ ) 

Nr - 

ip (da) (da) . 



2167r 
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This completes the decomposition of S^zw{U,v,a}. The renormahzation group 
invariant part is the Wess-Zumino-Witten term relevant for the effective theory 
at fixed A^: 

Jdx C^^w ^ ^wzw{U,v,a} . (C--!) 

The extension to the degrees of freedom carried by the rj' contains the following 
additional contributions, which are gauge invariant: 

-iD^e {R^pDMU^ - UpD^UW) } . (C.2) 
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